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I. INTRODUCTION 


This Final Report summarizes the work performed under Contract No. 
NASw-1726 for the period 24 May 1968 through 24 November 1969. The over- 
all problems associated with the physics of planetary atmospheres and the 
broad variety of scientific areas which have been considered are discussed 
most conveniently under the following major categories: 

Experimental Studies 
Theoretical Studies 

During the performance of the subject contract , a number of reports 
were generated which have been submitted and/or accepted and/or published 
in accredited journals. In some cases, a more detailed discussion of the 
material is available in the form of published GCA Technical Reports. In 
any event, the pertinent information concerning the journal publication is 
presented below for each of the categories previously cited. 

A. Experimental Studies 

"Primary Processes In the Photolysis of SO- at 18498," P. W'araeck, 
et al. J. Chem. Phys. 72, 3736 (1968). 

"Photoelectron Spectroscopy of the Rare Gases," J. Samson, et al. 

Phys. Rev. 173, 80 (1968). 

"Scattering Cross Sections of Argon and Krypton in the VUV, " J. Samson, 
et als. J. Opt. Soc. Am. 58 (1968). 

"Photo-Attenuation Cross Sections of Xe and Xe 2 between 1050 and 155oS," 
Shardanand. J. Quant. Spectrosc. Radiat. Transfer 8, 1373 (1968). 

"Angular Distributions of Photoelectrons and Partial Photo ionization 
Cross Sections," J. Samson. Proc. Roy. Soc. (in press). 

"Absorption Cross Sections of Minor Constituents in Planetary Atmos- 
pheres from 1050 - 2100A,” J. A. R. Samson, GCA Technical Report No. TR-69-7-N. 

"Angular Distribution of Photoelectrons,” J. Samson. J. Opt. Soc. Am. 

59 (1969). 

"Conversion Table for Wavelengths to eVs and cm" 1 ," J. Samson. GCA 
Technical Report No. TR-69-2-N. 

"Formation of N?* in Nitrogen," P. Wameck. J. Geo. Res. 74, #1, 

1396 (1969). 4 

"Higher Ionization Potentials of Nitric Oxide," J. A. R. Samson. Phys. 
Letters, 28A, #5, 391 (1968). 
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"Simultaneous Photo excitation and Photoionization of Helium, " J. Samson 
Phys. Rev. Letters 22 , 693 (1969). 

’Vacuum Ultraviolet Absorption Cross Sections of CO, HG£, and ICN 
between 1050 and 2100$,” J. Myer and J. Samson. J. Chem. Phys. 52, 266 (1970). 

"Heats of Formation of CH0"*and C3H3 by Photoionization, " C. Matthews 
and P. Warneck. J. Chem. Phys. 5JL, 854 (1969). 

"Absorption Cross Section and Photo ionization Yield of between 1050 
and 2200S, " J. Myer and J. Samson. J. Chem. Phys. 52 , 716 (1970). 

"Photon Scattering by Argon in the VUV," R. Cairns, F. Marmo and J. 
Samson. J. Opt. Soc.,Am. 6_0, 211 (1970). 

"Electronic States of CgH^" (J* Samson) Chem. Phys. Letters 4, 247 (1969). 

"Line Broadening in Photoelectron Spectroscopy," Rev. Sci. Instr. 40, 

1174 (1969).' 

"On the Measurement of Rayleigh Scattering," J. Samson. J. Quant. Spec- 
troscopy and Radiative Transfer 9 _, 875 (1969). 


B. Theoretical Studies 


"The Effect of Oxygen Cooling on Iono, spheric Electron Temperatures," 

A.' Dalgarno, et al. Planet. & Space Sci. 16 , #11, 1371 (1968) 

"Dipole Properties of Molecular Hydrogen," G. A. Victor and A. Dalgarno 
J. -Chem. Phys. 50, #6, 2535 (1969). 

"Mariner 6: Origin of Mars Ionized Carbon Dioxide Ultraviolet Spectrum,” 
A. Dalgarno, T. Degges, and A. Stewart. Science 16 7, #3924, 1490 (1970). 

"Electron Impact Excitation of the Dayglow," A. Dalgarno. J. Atmos. 

Sci. 26, #4, 753 (1969). 

"CO2 Dayglow on Mars and Venus,’* A. Dalgarno and T. Degges. Proc. I.A.U. 
Symposium 40 on Planetary Atmospheres, Marfa, Texas (October 27, 1969). 

"Carbon Atoms in the Upper Atmosphere of Venus," F. F. Marmo and A. 
Engelman, ICARUS 12, 128-130 (1970). 
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II. SUMMARY. OF TECHNICAL WORK PERFORMED FOR THE 
PERIOD 24 MAY 1968 THROUGH 24 NOVEMBER 1969 

The technical progress accomplished under the program is described 
in the present section in two general categories: (A) experimental studies 
and (B) theoretical studies. Each major category in turn is presented 
in two subsections, namely ( 1 ) published data, and ( 2 ) unpublished data. 
The first category contains all data which have been published in the open 
literature and are available to the interested reader. Accordingly, in 
these cases the actual published abstracts are repeated in the present 
section for general information and scientific content. The category of 
unpublished data contains descriptions of the studies performed which have 
not as yet appeared in the open literature. Accordingly, in these cases, 
it is considered appropriate to offer more detailed discussions since 
the data are not readily available elsewhere. 

A. Experimental Studies 


Published Data . 

"Heats of Formation of CHO and C^H^ by Photoionization" 

C. Matthews' and P. Warneck 
Journal of Chemical Physics, 51, #2, 854 (1969) 

• The ions CHO* and C 3 H 3 + occur with significant yields in the mass 
spectra of many oxygenated and unsaturated hydrocarbons, respectively, so 
that their heat of formation is of general interest. A particular need 
for the thermochemical properties has arisen from observations indicating 
that CHO + and C 3 HJ are major ions in the chemi-ionization of hydrocarbon 
flames (Refs 1-3). For a specification of the detailed reaction sequences 
generating these ions, reliable values for their heat of formation are 
required. 


"Formation of N^ in Nitrogen" 
P. Warneck 


Journal of Geophysical Research, Space Physics 
74, #1, 396 (1969) 


This paper is concerned with the discrepancies reported in 
ature over the preceding 15 years on formation in nitrogen, 
whether the process 


N 2 + N 2 ■* 


the liter- 
specifically. 


is a two body or a three body recombination process. Reports by Fite, et al . 
(Ref. 4) Saperoschenko (Ref. 5), Varney (Ref. 6 ), and Krovar, et al . (Ref. 7) 
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had indicated a two body process, whereas the more recent data by Warneck 
(Ref. 8) and Asundi, et al. (Ref. 9) displayed third order kinetics. The 
photoionization mass spectrometer technique established by Warneck (Ref. 10) 
should give straightforward results, but to counter a criticism voiced by 
Varnet (Ref. 11) new experiments were performed. The ionizing radiation 
used in these experiments was centered at 790&, close to the threshold of 
N 2 ionization, in order that the formation of vibrationally excited N2**" ions 
was precluded. Previous experiments had been performed with 764& radiation. 

Rate coefficients for formation were determined in one series of 
experiments in which the ratio of electric field strength to pressure, E/P, 
in the reaction region was held constant for several field strengths and 
pressures; and in another series of experiments, the pressure was varied 
while the field strength was held constant but differed for each set of 
values. All observations showed a pressure dependence of the derived bi~ 
molecular rate constant similar to that found previously (Ref. 10) thereby 
demonstrating third order kinetics. Moreover, the third order rate co- 
efficients derived from all these experiments agree within the experimental 
error, the value being k = (4.4 +0.6) x 10*29 cm^/molecule^ sec. Thus, 
the variation with pressure is due to the influence of .third body and not 
due to a variation in E/P in the region of field strengths covered (8 < 

E/P < 33 volts/cm torr>. 

"On the Measurement of Rayleigh Scattering” 

J.A.R. Samson 

J. Quant. Spectrosc. Radiat. Transfer £, 875(1969) 

The general expression for the angular distribution of Rayleigh scattered 
radiant energy from a gas is developed for the case of partially plane 
polarized incident radiant energy. It is shown that, if measurements are 
made at specific angles, then both the total scattering cross section and 
the normal depolarization factor can be determined without knowing the de- 
gree of polarization of the incident radiant energy. 

"Line Broadening in Photoelectron Spectroscopy" 

J.A.R. Samson 

Review of Scientific Instruments 40, #9, 1174 (1969) 

The factors causing line broadening in photoelectron spectroscopy are 
discussed. Experimental results are given for the half -width of the 584? 
helium line emitted from a variety of sources. Typical half-width values 
lie between 1 and 6.5 mV, which corresponds to 0.027 and 0. 18.S, respectively. 



Curves are presented showing the contribution of the thermal motion of 
the photoionized gas to the spread in- the energy of the ejected photo- 
electrons. 


"Electronic States of C^Hg’ 1 


J.A.R. Samson 

Chemical Phys. Letters 4, #5, 257 (1969) 

The photoelectron energy spectrum of benzene has been measured with 
monochromatic radiation of wavelengths between 1253 and 584A. The first 
two ionization potentials are shown to be 9.25 and 11.5 eV. Explanations 
are given which account for previously observed structure lying between 
the first two ionization potentials. 

"Simultaneous Photoexcitation and Photoionization of Helium 
J.A.R. Samson 

Physical Review Letters 22, #14, 693(1969) 

The probability of simultaneously photoexciting and photoionizing 
helium has been measured near threshold at a wavelength of 186A. It was 
found that about 8+2 percent of the ions formed at that wavelength were 
excited into the n = 2 state. 

"Angular Distribution of Photoelectrons" 

J. A. R. Samson 

Journal of the Optical Society of America 59, #3, 356 (1969) 

Transition probabilities for transitions from the ground state to the 
ionic state of atoms or molecules can be determined by using the technique 
of photoelectron spectroscopy. This technique measures the kinetic energy 
of the ejected photoelectrons and the relative numbers of these electrons 
within each energy group. Knowing the relative numbers of electrons in each 
group and the total photoionization cross section at a specific wavelength 
we can obtain the transition probabilities at that wavelength. However, 
it is important to obtain the true representation of the relative numbers 
of electrons within each energy group. 

"Higher Ionization Potentials of 
Nitric Oxide" 

J.A.R. Samson 

Physics Letters 28A ,#5, 391 (1968) 
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Higher ionization potentials of nitric oxide are reported using the 
technique of photoelectron spectroscopy. Previously unreported ionization 
potentials of 21.72 and 23.1 eV were observed. 

'Vacuum Ultraviolet Absorption Cross Sections of CO, HC0 , and ICN Between 

1050 and 21008" 

J.A.R. Samson 

Journal of Chemical Physics 52 , #1, 266 (1970) 

Absorption cross sections have been measured under intermediate re- 
solution (0.258) for CO, H C& , and ICN in the spectral region from 1050 to 
2100S. An attempt has been made to analyze the spectra in terms of the 
reactions responsible for the observed features* A Rydberg series has 
been identified in the absorption spectrum of ICN with a series limit of 
1140.6 ± 0. l8. 


"Photoelectron Spectroscopy of the Rare Gases" 

J.A.R. Samson 

Physical Review 173 , #1, 80 (1968) 

A spherical retarding potential analyzer was used to measure the 
kinetic energies of electrons ejected by photoionization of the rare gases 
as a function of wavelength from the ^q/2 ionization threshold down to 
about 4008. The ratio of the number of ions produced in the ground z P 3 / 2 
state to the number produced in the excited s ^ a ^ e was f° un d to be 

constant with respect to wavelength and had the following values: 2.18, 

1.98, 1.79, and 1.60 for Ne, Ar, Kr, and Xe, respectively. From these 
data the specific photo ionization cross sections for the rare gases were 
obtained. At wavelengths shorter than the MU, Nj_, and 0-^ thresholds for 
Ar, Kr, and Xe, respectively, the ejection or s-electrons was observed 
•and -the specific cross sections for their production measured. 

"Absorption Cross Section and Photoionization Yield of I2 Between 1050 6e 2200i 

J.A.R. Samson 

Journal of Chemical Physics 52 , #2, 716 (1970) 

The total absorption cross- section and photoionization yield of mol- 
ecular iodine were measured from 1050 to 220o8. Two thresholds were ob- 
served in the photo ionization yield curve. The first was- associated with 
the ion-pair production 1*^ + I” and the second to direct ionization of 12* 
From the appearance potential for ion -pair production and the accurate 
ionization potential of I 2 it was possible to calculate the value' of the 
electron affinity of iodine, namely, 3.073 + 0.014 eV. 
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"Primary Processes in the Photolysis of S0 2 at 1849$" 

J-.N. Driscoll and P. Warneck 
Journal of Physical Chemistry 12, 3736 (1968) 

The production of SOg in the 1849$ photolysis of pure S0 2 and in mix- 
tures containing 5-75 percent S0 2 in 0 2 was studied at atmospheric pressure. 

An interference filter isolated the 1849$-1924$ lines from a mercury arc 
and eliminated the emission at longer wavelengths. The quantum yield found 
for SOg formation in pure S0 2 was 0.50 + 0.07. The addition of oxygen 
increased the SOg quantum yield to a plateau value of 1.04 + 0.13, but 
for mixtures containing less than 30' percent SO 2 the quantum yield declined 
again. The photodecomposition of S0 2 at' 1849$ is complicated by the occurrence 
of several primary processes including direct dissociation due to an absorp- 
tion continuum, excited S0 2 molecule formation, and predissociation. The 
observed SOg quantum yield was utilized to evaluate the relative importance 
of the various primary processes. 

"Photon Scattering by Argon in the Vacuum Ultraviolet" 

R. Cairns, F. Marmo & J. Samson 
Journal of the Optical Society of America 6j0, #2, 211 (1970) 

Scattering cross sections were measured with the method of Gill and 
Heddle, but with dispersed incident flux. These measurements constitute 
a sensitive test for the presence of rare-gas quasimolecules. The ratio 
of At 2 to Ar at 100 torr was deduced to be 8.6 x 10“5. At 1090$ the 
molecular cross section is estimated to be 1,6 x 10~^-^ cm^. The attenua- 
tion cross section is 8.0 x 10"^*® <^2, 

Unpublished Data . -"High Resolution. Quantitative Absorption Data of Minor 
Constituents in Planetary Atmospheres." The fact that high-resolution data 
are of significant importance in isolating the contribution due to minor 
constituents in a planetary atmosphere has been amply demonstrated in the 
infrared spectra of planetary atmospheres. For certain specific cases, 
similar potential exists in the VUV region so long as the appropriate data 
are available. For example, it is propitious to perform measurements in 
those vacuum ultraviolet (VUV) spectral regions where the minor atmospheric 
constituent displays a cross section which far exceeds all others assoc- 
iated with the major constituents. These features can be even more easily 
isolated if the absorption features are discrete. Performance of this 
laboratory task is facilitated by employing a new special McPherson vacuum 
monochromator which GCA has recently acquired. A description of the salient 
features and the performance capabilities of this monochromator is included 
in the discussion below. 



During the current period , the general experimental capability of 
obtaining VUV high resolution spectra was investigated. For these pre- 
liminary measurements, two gases were examined, namely methane and benzene. 
The results of these investigations are given below. 

The monochromator is a McPherson model No. 225 which is a one -meter 
normal incidence instrument equipped with a 1200-^ /mm grating. The mono- 
chromator is capable of resolving 0. 2& in the first order. Attached to 
the exit slit of the monochromator is a double.-beam system which consists 
of a wedge shaped mirror oscillating a 6 cps. The radiations from the 
monochromator strikes each face of the wedge at a grazing angle and is 
reflected into -a reference and a sample cell. The surfaces of the wedge 
are concave which focuses the radiation into the absorption cells and 
prevents loss of light due to a diverging beam (see Figure 1). The cells 
are completely enclosed with LiF windows in front and in back of the cells. 
After passing through the cells, the radiation impinges onto sodium sal- 
icylate screens converting the vacuum UV radiation into near UV radiation 
centered around 4100&. This radiation is then detected by two matched 
photomultipliers. The most important and difficult problem is the match- 
ing up of the two pairs of LiF windows in both the reference and sample 
cells. Once this has been achieved, the log-ratio of the output signals 
of the two photomultipliers is proportional to the total absorption cross 
section of the gas in the sample cell. 

The McPherson log-ratiometer is skillfully designed to allow log- 
ratios to be made with the many- lined spectrum of H2* In order to facil- 
itate the ratioing of rapidly varying intensities over several orders of 
magnitude in intensity, the unit first takes the signal i r and if from 
the reference and sample-cell photomultipliers, respectively, and feeds 
the signals into a time sharing log amplifier. The outputs of the log 
amplifier (log i r and log if) are then fed into a difference amplifier 
giving an output 


log i r - log if = log (i r /if) 

That is, the output gives the logarithm- to the base ten of the ratio 
(i r /if), (See Figure 2) which is related to the tota-l absorption cross 
section as follows: 


1,123 log (ir / i£ ), 

where n is the number density of the gas in the sample cell and L is the 
length of the sample cell. 

Preliminary results: 

The total absorption cross sections of methane and benzene have been 
measured from 1100 to 2100$. The results are shown in Figures 3, 4, & 5. 
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Methane. . 


i 

7The only previous cross-sectional data which exists in -the spectral ■ 
region for methane is that'' of Sun and Weissler (Ref. 12) and of Watanabe, 
et ai . (Ref. 13). Weissleiis data has a long-wavelength limits of 1300$ 
and ^furthermore was taken at only a. few. points with a line spectrum. Thus, 
the data of Watanabe, et al. in -1953 represents the only comparable data 
with the presently reported results. However, their wavelength resolution 
was,, only 1$ whereas the present data was obtained with a resolutibn of 
0.45$., Watanabe 's -absolute cross sections agree with the present’ results 
over :most of the wavelength region. However, the present data resolves . 
structure below 1125$, which they could not observe. With the exception- 
of this structure, the absorption of methane is continuous from d?out 
1450$ to shorter wavelengths. No absorption was observed between 1450$ 
and 2100$. 

Benzene 

The benzene spectrum, nas oeen scuaiea tnorougniy „-Dy wiiKinson liter, rqj 
using photographic, techniques at .high, resolution.. However, total absorptior 
cross sections have only been made in a few, isolated wavelength regions 
or else made at low wavelength resolutions The present data shown in_„ 
Figures 4 and 5 shows considerably, detailed absorption structure not read- 
ily seen with poor resolution ..-'•"'In Figure 4, two Rydberg series, denoted 
at R and R 1 , can be seen leading up to the first ionization potential qf 
benzene at 1340$ ^ ' ' 4 

"Drift Velocity -Measurements for Atmospheric Ions.-" In the course of 
investigating drift yelocities of N 2 ^ ions formed by photoionization of * 
’nitrogen, it was found that when the ionizing wavelength was 790$ ? , the 
observed drift velocities were in excellent agreement with data obtained 
by other experimenters .using techniques different from ours; however, when 
the spectral setting -of the monochromator was 764?, the drift velocities 1 
observed were by about* a factor of two greater. Simultaneously, tiie ion 
density profile observed on the screen of the oscilloscope was noj -longer 
symmetric, but became skewed, indicating the presence of two groubs of 
ions 'with different Velocities. At the wavelength setting of 790$ (close 

to the threshold of N 2 + "~fprmation) only ground- state N 2 * ions can, be gen- 
erated. At 764$ excitat iofKto the first vibrational level is also possible. 
This! channel would be significant if an appreciable fraction of the ions 
were. formed by autoionizatioji. The observation that two groups of ions 
are involved at 764$ makes" evident that autoionizatioh * is significant and 
that,.N 2 *^ (v = l)_JLons~-are endowed with a higher mobility than- N^lCv = 0) 
ions. 

fHe ~ f orma t idnlof occi^T“a^^dtE“ wavelength s ! ^'THe'ra t e " c ; o ef f ici en t 
for this reaction had been measured previously with the photoionization 
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technique at 764$. Third body kinetics were observed. This observation 
has been criticized as being inconsistent with a number of other data in 
the literature. .In view of the participation of N 2 + (v = 1 ) ions at 764$ 
it was considered necessary to re-investigate the process of N 4 + forma- 
tion at 790$ where only ground state N 2 * ions are initially present. The 
details of the results and their implications are given in the present 
section in the form of a paper submitted for publication in the Journal of 
Geophysical Research. - 

In a recent comment concerning a discussion (Ref. 8 ) of atmosphere 
ion neutral reaction rates, Varney (Ref. 6 ) objected to listing the forma- 
tion of via the association 

4 + N 2 - N+ (1) 

as a three body reaction, and he cited- experimental results which indicated 
that it proceeds as a two-body process. Specifically , the after glow data 
of Fite, et al . (Ref. 4), the mass spectrometer results of Saporoschenko 
(Ref. 5), and variety of drift experiments (Refs. 6 , 7, 15) have indicated 
a direct two -body association whereas the more recent photoionization mass 
spectrometer experiments by Warneck (Ref. 10) and the low pressure electron 
impact mass spectrometer observations by Asundi, Schulz, and Chantry (Ref. 

9) displayed third order kinetics. Since the discrepancy is serious, we 
have performed a new set of experiments with the photoionization technique, 
following the suggestions by Varney (Ref. 11) that the parameter E/p, the 
ratio of the electric field strength to the pressure, be kept constant in 
the ion source so that the energy of the reacting ions remains unchanged. 

The new results and their implications are the topic of this note. 

Initially, the former results with constant E and varying p were ver- 
ified, using for the production of N 2 + radiation centered at 764$, as des- 
cribed previously (Ref. 10). Since this permits the production of vibra- 
tionally excited nitrogen ions, the new experiments were performed at 790 $, 
near the threshold of nitrogen ionization, where the formation of excited 
nitrogen ions is precluded. At 790$ the measured drift velocities of N 2 + 
agreed well with those reported in the literature (Refs. 15, 16) whereas 
the drift velocities at 764$ were by a factor of two higher, presumably 
-due to the presence of a considerable portion of vibrationally excited nitro- 
gen ions formed by autoionization. Rate coefficients for N 4 * formation were 
determined at 790A in a series of experiments for which E/p was held con- 
stant for various field strengths and pressures; and in another series of 
experiments for which the pressure was varied while the field strength 
was held constant but differed for each set of values. All observations 
showed a pressure dependence of the derived bimolecular rate constant sim- 
ilar to that found previously (Ref. 10), thereby demonstrating third order 
kinetics. Moreover, the third order rate coefficients derived from all 
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these experiments agree within the experimental error,, the value being 
k = (4.4 + 0.6) x 10“29 cc^/molecule^ sec. Accordingly, the variation with 
pressure is due to the influence of the third body and not due to a variatior 
with E/p in the region of relative field strengths covered ( 8 < E/p < 33 
volt/cm/torr) . The new results, therefore, corroborate our previous con- 
clusion that for the conditions employed in the photoionization mass spec- 
trometer an equilibrium between N 2 *** and is not established. 

The ionizing energies employed here in photoionizing nitrogen are in a 
region where corresponding electron impact experiments have demonstrated 
formation from neutral excited nitrogen via the process: 

e + N 2 N* + e 
N* + N 2 -> N* + e 

occurring in addition to N 4 "** formation from N 2 + (Refs. 9, 17) 

In view of the knowledge that in the wavelength region 658$ - 850$ the 
absorption spectrum of nitrogen, is at least partially discrete, and that 
photo ionization efficiencies are less than unity, the formation of excited 
nitrogen molecules is certain to occur also in the present experiments. 
Although the excited states reached by UV radiation and by electron impact 
are not necessarily the same, it is clear that allowance must be made for 
the possibility of reaction (3) participating in the photoionization mass 
spectrometer experiments. However, no experimental evidence has been found 
for any occurrence of the process. Specifically; no N 4 ~ f ' formation was de- 
tected at energies below (i.e. at wavelengths above) the onset of the N 2 * 
formation — in contrast to the electron impact results, and approximately 
the same rate of N 4 * formation as at 790$ was observed when the wavelength 
wa,s set to 685$, where the ionization efficiency is nearly unity. Thus, 
while reaction ( 3 ) cannot be entirely precluded, it does not appear to 
occur to a significant extent in the present experiments. It is concluded 
that the observed formation of N^**" is due predominantly to reaction ( 1 ) in 
its termolecular mode, in agreement with our earlier conclusion (Ref. 10). 

There remains to explain the discrepancy arising from the various data 
in the respect to third or second order kinetics of N 4 + formation. Evi- 
dently, this problem is connected with the observation that at high E/p 
values ions gain sufficient energy in the field to dissociate upon col- 
lision. Mass spectrometrically identified ion transients in drift tubes 
are nearly identical for N 2 * and when E/p is high, demonstrating that a 
rapid equilibrium between the two species is achieved (Refs. 15, 16). In 
low feilds, however, the two ion transients are well separated even though 
a significant N 2 + - N 4 + conversion occurs. In this region, therefore, a 


( 2 ) 

(3) 
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rapid interchange is not operative. The low field region is well suited 
to studies of reaction (1) and in this region., the formation of is . 
third body dependent. On the other hand, in the high field region, reaction 
(1) is in equilibrium with the dissociation of so that in this case, 

the abundance of will depend only linearly on the nitrogen pressure. 

It appears that in those cases where bimolecular N, + formation was deduced 
(Refs. 4,5,15) the high field conditions favored 'equilibration of N 2 + 
and N 4 + . The mass spectrometer results of Asundi, Schulz and Chantry (Ref. 9) 
provide direct evidence for such an interpretation. 

"VUV (X\ 1050-2000£) High Resolution Absorption Cross Section Data with 
Emphasis on Minor Constituents In Planetary Atmospheres.” With reference to 
the statement of work under this contract, Task Item #3 has been completed. 
The VUV absorption cross sections have been measured for a number of gases 
that represent possible atmospheric constituents. In most cases high resol- 
ution absorption spectra have not been heretofor obtained in this spectral 
range, in spite of the fact that such data are directly applicable to the 
determination of planetary atmospheric constituency. For example, such 
data could be of particular value in those spectral regions where the ab- 
sorption properties of the major constituents are relatively weak and/or 
constant whereas the minor constituents display strong discrete absorption 
features due to resonance bands, Rydberg series, etc. In addition, high 
resolution VUV absorption data are important in the interpretation of the 
physical mechanisms responsible for the absorption processes. For example, 
such data provide information which can define the produce radical of a 
given photodissociation, the energy distribution of the products, the possi- 
bilities of their reacting with other species to form new products, etc. 

These investigations were accomplished with a McPherson Model 225 one-meter 
VUV scanning monochromator equipped with a McPherson 630 VUV light source. 

The optical features are described elsewhere (Ref. 18) so. that they are not 
repeated here. However, the unique feature of the experimental arrange- 
ment is the employment of a McPherson Model 665 double beam attachment 
that was incorporated with the monochromator. 

The unique feature of the present experimental arrangenent is a 
McPherson Model 665 double beam attachment that was incorporated with the 
monochromator. This attachment has two optical paths with identical absorp- 
tion cells (one is a reference cell) and a vibrating mirror that alternately 
reflects the dispersed light beam through each cell. The alternating sig- 
nals are measured by two independent photomultiplier tubes (EMI 9635 B 
With sodium salycylate coatings) and the logarithm of the ratio of the two 
signals as evaluated by a logarithmic ratiometer and then registered on a 
linear strip chart recorder. If a known gas pressure is introduced into 
one of the absorption cells whose length is accurately known, then the ab- 
sorption cross section in units of cm^ can be determined directly from the 
recording. 



For this technique to be most effective, it is necessary for the optical 
paths including the window transmissions of the cells to be as identical as 
possible at all wavelengths considered. When the optical paths are properly 
. balanced then one will obtain a smooth horizontal base line for effective 
zero cross section when there is no absorbing gas in either cell. 

The gases and vapors to be studied were introduced and controlled with 
a gas filling system which included a needle leak valve, a diffusion pump and 
several vapor traps, a McLeod gauge, a thermocouple gauge (0-100 microns), 
and an oil monometer for pressures up to 30 torr. Vapors were introduced 
through a double stopcock-flask arrangement on the same system. The gas 
pressures for the cross section measurements were read with the oil mono- 
meter over the range from 0.05 torr - 30.00 torr (+ 0.05 torr) which was 
adequate for most of the measurements. 

t 

The gases investigated included CO, NH 3 , CH^, NO, S0 2 , H 2 0, C0 2 , 0 2 , HC0 , 
12 and ICN; the resultant data are shown in Figures 6-19 inclusive. 

The data on the several gases shown in the figures have been discussed, 
analyzed, and interpreted in terms of the applicability to atmospheric 
physics in a GCA technical report, i.e. GCA TR 69-7-N by James A. R. Samson 
and James Myer. The interested reader can consult the cited reference so 
that it is unnecessary to include any further discussion on these investiga- 
tions . 

"Ion Molecule Reactions in Hydrocarbon Systems . 11 From spectroscopic 
observations, it is known that the Jovian planets contain appreciable amounts 
of methane and ammonia in their atmospheres, in addition to hydrogen and 
helium. McNesby (Ref. 19) recently has studied the photochemistry of such 
a system and has shown that, with hydrogen in excess, the first order photo- 
dissociation products of methane and ammonia essentially recombine, result- 
ing in no net change of the system except to heat the atmosphere. Neverthe- 
less, a small amount of higher hydrocarbons would be formed as a second 
order perturbation. For example, Sagan (Ref. 20) has argued for the forma- 
tion of adenine and other amino products. With this background, 'it becomes 
of interest to investigate ion reactions involving hydrocarbons. A great 
number of such reactions has been studied in the past with conventional 
■ electron impact mass spectrometers. With the applied techniques, these 
investigations have been complicated by the many ions produced simultaneously 
in the ion source, and their reactions with the parent gas. It appears, 
that the photoionization mass spectrometer technique, by virtue of the sel- 
ective ionization processes involved, • can reduce the complexities of the 
reaction system sufficiently to allow a less ambiguous analysis of the data. 

On this basis, we have studied ion reactions in the systems ethylene- 
acetylene, acetylene-methane and methane -ethylene. The results have already 
lead to some interesting conclusions, although all the data have not yet 


18 



ABSORPTION CROSS srcriON (Mb) 



ABSORPTION CROSS SECTON(Mb) 



c«o<5 sttnc/l {MbJ 



MGT*' [*) 


Figure 7. Absorption cross section of NH3 as a function of wavelength between 2200R and 1050S. 
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Figure 14 Absorption cross section of 0? as a function of wave- 
■ ■ ‘ length between 2000 and 1750 K. 
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Figure 15. Absorption cross section of HCU as a function of wavelength between 2100X and 1050X 
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Figure 17 . Absorption cross section of I2 as a function of wave 


length between 2020& and 1760&. 
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Figure 19. Absorption cross section of ICN as a function of wavelength between 3000 and 2100X. 
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been analyzed, and a detailed evaluation of those data which have been 
analyzed has not yet been possible. The results presently available, and 
their implications will be discussed below. 

The techniques employed in this study have been described in detail 
previously (Ref. 21). briefly, the apparatus involved a 180 degree magnetic 
mass spectrometer with conventional ion source geometry and electronic multi- 
plier. Ions are produced by extreme UV radiation generated in a repetitive 
spark light source. A Seya monochromator provides wavelength isolation with 
a spectral resolution of about 5&. The light beam transversed the ion chambe 
midway between the repeller plate and the ion exit orifice which were placed 
0.6 cm apart.' The repeller potential was one volt. The gas pressure in 
the ion source was measured with a McLeod gauge through the hollow stem of 
the repeller. 

Reactions were investigated in the low pressure domain where the motion 
of ions is mainly acceleration in the repeller field. Products were assessed 
from the appearance of new peaks as the pressure in the ion chamber was 
raised. Rate coefficients were determined from the decrease of the primary 
ion intensity as a function of the gas' pressure, or as a function of the 
partial pressure of an added second gaseous component; and from the re- 
sidence time of the primary ions in the source chamber. Residence times 
were determined from the time delay of the ion pulse arriving at the mass 
spectrometer collector after the deposition of ions in the source by the 
light pulse. As described previously, delay times were measured for sev- 
eral repeller potentials and a graphical extrapolation toward infinite 
repeller potential is used to separate the ion source residence time from 
the total flight time of the ions in the mass spectrometer. In the present 
experiments the ha If -rise time of the integrated ion pulse was utilized 
for reference in the delay time measurements so that averaged residence times 
were obtained. 

In the following, results will be presented for ion-molecule reaction 
in ethylene, acetylene, and their mixtures. At ionizing wavelengths above 
945A up to the ionization threshold, the primary photoionization products 
in both ethylene and acetylene are the parent ions C2H4 + and C 2 H •£ , 
respectively, and other ions are formed only by secondary reactions of these 
ionic species. In the wavelength region 900-945^, photoionization of acet- 
yl en ^ still yields only the parent ion, whereas that of ethylene produces 
C 2 H 2 in addition to C 2 H^*^. This latter wavelength region thus facilitates 
a study of the reactions of C 2 H 2 + with ethylene in the absence of inter- 
ference by acetylene. The present * investigation made use of radiation at 
1087A and 923A both provided by operating the spark source with nitrogen. 

The development of the ion spectrum with pressure when acetylene is 
ionized at 1087A is shown in Figure ,20. The important product ions appear 
at mass numbers 50> 51, and 52. The reaction paths responsible for the 
production of these ions are 
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c 2 H+ + c 2 h 2 


C 4 H 2 

C 4 H 3 " 

c 4 « 4 4 


+ products 
+ products 


(*> 


At pressures below 10 microns, the product ratios are found to be constant. 

The determination of the average values for the branching ratios of reaction 
(4) gave 29 percent for the formation of C4H2+ , 66 percent for C4H3 + and 5 
percent for C4H4+. The pressure dependence for the formation of the last 
ion is such that it precludes the participation of the third collision partner. 
The formation of 04X14+ in a binary collision is possible only if the life- 
time of the intermediate collision complex is sufficiently long, so that 
the collision complex is stabilized by subsequent collisions regardless of 
the pressure. Since the produce branching ratio does not depend on the 
pressure, it must further be concluded that the other two products result 
from the same process as C4H4 + . Thus, the detailed mechanism of reaction 
(4) in the investigated pressure domain can be written as 


C 2 H 2 +C 2 H 2 


C 4 H 4 + * + C 2 H 2 


+* 


C.iL 
4 4 


C 4«2 + + C 2»4* 

c 4 h 3 + + c 2 h 3 * 


(4a) 

fcb) 


C,H, 
4 4 


where the asterisk denotes (vibrationally) excited species. The present- 
results contain no information on the fate of the neutral products in 
reaction (4b)., but a possible course is their dissociation to yield H2 or 
H and acetylene. 


At pressures greater than 10 microns, the ratio of products H50/M51 
decreases and new ions appear at mass numbers 74, 75 and 76 due to the 
reactions 


C 4 H 2 +C 2 H 2 


C 4«3 +C 2 H 2 


c 6 h 2 + + a, 
C 6 H 3 + + H 

C 6 H 3 + + H 2 
C 6 h 4 + + H 


(5) 


( 6 ) 
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of these , r eat ion (5) predominates, producing and CgH^ ions in approx- 

imately equal amounts. The sum of all 'the observed ion intensities, i.e. 
that of the acetylene and that of the various product ions, was found 
directly proportional to' the acetylene pressure, indicating that reactions 
(4) through (7) account for essentially all the ions observed and further 
reactions are negligible. The rate coefficient for reaction ^4) including 
all its channels can be derived from the decrease of the 0^2 i° n i- nten ~ 
sity as a function of pressure from the relation 


2.3 


N 

i°s sr 


where i is the average residence time of C2H2 in the ion source, n the 
number density of acetylene in the source, N Q is the sum of all ions taken 
equal to the initial C2H2 + ion intensity and N is the measured C2H2+ ion 
intensity. The average value of the rate coefficient thus determined for 
reaction (4) is k]_ = 1.39 x 10~ 9 cc/raolecule sec. The value expected on 
the basis of the ion induced dipole interaction theory is k]_ (theor.) = 

1.29 x 1CT 9 cc/molecule sec. The experimental rate coefficient is found 
to agree with the theoretical one within experimental error. 

Similar results were obtained for ion reactions in acetylene when the 
photoionizing wavelength was 923.S. Again, C2H2"*" was the only primary ion 
observed, and reaction (4) with its three channels was the major loss pro- 
cess for C2H2 + * The ratio of product channels was determined as C4H2+ 29 
percent, 04113+ 67 percent, C4H4* 4 percent. Within experimental error, these 
values are the same as those given above for the reaction when working with 
1087$ radiation. Moreover, the associated rate coefficient was found to 
be k x = 1.38 x 10“ 9 cc/molecule sec. This value agrees well with the value 
given above, and like that one is in accord with theoretical prediction. 

These results, therefore, show that the acetylene ion formed at 923A behaves 
identically to that formed at 108 7$, even though it may be produced with 
some internal excitation. 


The development of the spectrum- of the major ions with pressure when 
ethylene is ionized at 1087$ is shown in Figure 21. The significant pro- 
duct ions appear at mass numbers 41 and 55. They are due to the reaction 


C 2 H 4 + + C 2 H 4 4 C 3 H 5 + + CH 3- (7) 

- C^ + H 

The branching ratio of this reaction is found to result to 90 percent 
in C3H5+ ions, and to 10 percent in C4H7+ ions. At pressures above about 
10 microns, this ratio is perturbed by the strong increase of a tertiary pro- 
duct at mass 69 due to the condensation reaction 
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ETHYLENE PRESSURE (microns ) 

Figure 21 Intensities of ^jor ions in ethylene as a function of pressure, 
lonxzing wavelength 108 7A. 




( 8 ) 


C 3 H 5 + + C 2 H 4 - C 5 H 9 + 

which consumes the ions formed in Reaction (4). The C 4 H 7 ’*" ion formed 

in the reaction (Ref. 22) also undergoes further reactions but these are 
significant only at pressures greater than 30 microns and even then their 
contribution to the total ion abundance is small. Contributions to the 
total ion spectrum were observed also on mass numbers 42 and 70, but they 
were too small to be taken into account in a meaningful way. However, 
they are observed with higher intensity in the photoionization of ethylene 
at 923& and will be discussed further below. 

As in acetylene, an association product, C^Hg***, was found in the 
photoionization of ethylene. However, compared to the behavior at the mass 
54 ion in acetylene, .the association product in ethylene showed a differ- 
ent pressure dependence. Figure 22 indicates that the ratio of the inten- 
sity of the mass 56 ion to the total product ion intensity is not constant 
but increases with pressure. The dependence on pressure is linear and 
the straight line obtained extrapolates to the origin of the plot. This 
is evidence that in contrast to the acetylene case, the association product 
of the Chylene ion is formed by a third body reaction 

C 2 H 4 + + 2 C 2 H 4 •* C 4 H 8 + + C 2 H 4 (9) 

The rate coefficient for reaction (7) was determined from the decay of the 
C 2 H 4 * ion intensity compared to the total ion intensity. A value of = 
1.24 x 10"9 cc/molecule sec was obtained. The rate coefficient calculated 
from the induced dipole measurement theory is k^Ctheor. )=1. 29xl0 - 9 . As 
for the previous cases, the agreement between the experimental and the 
theoretical values is good. From the rate coefficient associated with 
reaction (7) and the data shown in Figure 22, the rate coefficient for 
reaction (9) can be derived as kg = 1.5 x 10"25 cc^/molecule^ sec. The 
magnitude of this value indicates that the average lifetime of the inter- 
mediate collision complex before it is stabilized by a third body collision 
is approximately 1 x lO”? seconds. Since this rate coefficient is derived 
from a product intensity rather than the primary ion intensity decrease, 
is is probably not as reliable as the other rate coefficients. Neverthe- 
less, the implied order of magnitude makes reaction (9) one of the most 
rapid termolecular ion reactions. Gordon and Ausloos, (Ref. 22) from 
studies of the photolysis of ethylene at 1048^/66^ had concluded that for- 
mation of the butenes is due mainly to neutralization of the C^Hg** - ion 
produced by reaction (9). 

Since the primary ions formed by photoionization of acetylene and 
ethylene with 1087.& radiation are only the parent ions of these gases, and 
their reactions have been established by the results reported above, it 
was of interest to look for additional reactions of these species in a 
mixture of ethylene and acetylene. Specifically, charge transfer from the 
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acetylene ion to ethylene was expected to occur, as it is an exothermic 
process. In these experiments, therefore-, acetylene was entered into 
the ion source chamber at a fixed rate providing a partial pressure of 
3.8 microns, and ethylene was added in varying amounts at partial press- 
ures up to 8 microns. In this pressure range, the ion intensities due to 
ethylene and its reaction products at mass number 28, 41 and 55 fitted well 
with the intensity curves shown in Figure 21 for pure ethylene so that the 
ethylene ions must react little or not at all with acetylene. The acety- 
lene ion intensities, however, decreased considerably with increasing 
ethylene pressure as shown in Figure 23a indicating the occurrence of an 
interaction^' At the same time new peaks appeared at mass numbers 39 and 
53. Thus the reaction 

C 2 H 2 + + C 2 H 4 


is indicated. From the intensity distribution of the ions M39 and M53 the 
branching probability for reaction (10) was obtained as 03113“*" 71 percent 
and’ C4H5* 29 percent. The sum of the intensities of these products matched 
the losses of the acetylene ions within the experimental error (see Figure 
23a) so that reaction (10) must be the major reaction that acetylene ions 
undergo with ethylene. Several previous workers (Ref. 23) have concluded, 
however, that charge transfer is a major reaction channel. To investigate 
this problem further, the sum of the intensities of ions due to C2H4'*' and 
its interaction with ethylene are plotted in Figure 23b for ethylene in 
the presence of acetylene and for ethylene along. It is evident that both 
sets of data coincide, whereas if charge transfer took place to an appre- 
ciable extent, the total ethylene ion intensity should increase with in- 
creasing ethylene pressure. It may be agreed that the lack of an increase 
in the ethylene ion intensity is due to a reaction of these ions with acety- 
lene and that charge transfer is operative nonetheless. The probability 
for this additional reaction would be greatest at low ethylene concentra- 
tions because then the concentration ratio of acetylene to ethylene is large. 
However, with increasing ethylene pressure the importance of this process 
diminishes in comparison with the interaction with ethylene, and the total 
ethylene ion intensity should increase correspondingly. From the failure 
to detect such an increase at the highest ethylene pressure involved 
(C2H4/C2H2 =2), from the intensity losses of the acetylene ions, and from 
the corresponding increase of the M3 9 and M53 ion intensities shown in 
Figure 23a it is estimated that the charge transfer reaction 

C 2 H 2 +C 2 H 4 - C 2 H 4 + C 2 H 2 (11) 

oc r,tT ' fi Hian 15 percent of the overall reaction of acetylene ions. 


- c 4 h 5 + +h 


( 10 ) 
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From the decrease of the 02^2^ i° n intensity as the ethylene pressure 
is increased and from the measured residence time of the acetylene ions in 
the source chamber, the rate coefficient for the reaction(lO) and (11) combir 
is found to be k 7 -fg = 1.35 x 10 “9 cc/molecule sec. The rate coefficient 
calculated from theory is k 7 -f$ (theor. ) = 1.31 x 10”9 cc/molecule sec. With! 
the experimental error the values are identical. 

When ethylene was photoionized with 923$ radiation, the major ion 
reactions were found to be reactions (7) through (10) discussed above, in- 
volving the primary ions C 2 H 4 + and C 2 H 2 + ‘ In accord with expectation the 
major product ions appeared at mass numbers 41, 55, 56, and 69, and at 
39 and 53 respectively. However, some additional reactions occurred also, 
as evidenced by signals on mass numbers 27, 29, 40, 42, 54, 67 and 70. Al- 
though the contribution of these ions amounted to no more than 7 percent of 
the total ion intensity, their assignment is required for a complete inter- 
pretation of the results. The M67 ion is taken to be the condensation pro- 

duct of C 3 H 3 ' } * with ethylene, because it was present also in the 1087^ 
photoionization of an acetylene-ethylene mixture when the ethylene pressure 
exceeded 10 microns. The behavior of most of the other ions identifies 
them as product ions. By analogy to reactions (7), ( 8 ) and (10), therefore, 
the signals on mass numbers 40 and 54 are assigned to products of M27 , and 

the signals on mass numbers 42 and 70 are assigned to products of the ion . 

M29. These assignments then reduce the problem of the minor ions to the 
origin of M27 and M29. The intensity of these ions and their follow-on 
products is much greater than the expected isotope contributions in ethylene 
(mainly due to C^*), so that the M27 and M29 ions must be identified as 
^ 2 ^ 3 * and C 2 H 5 **'. They give rise to the reactions 


C 2 H 3 + C 2 H 4 


C,H* +. CH 0 
3 4 3 


. C.H t + H 
4 6 


( 12 ) 


From the M40 and M54 contributions the branching ratio of this reaction is 
found to yield 03 ^^ to about 60 percent and to about 40 percent. Sim- 

ilarly, M42 and M70 are due to 


C 2 H 5 +C 2 H 4 
C 3 H 6 + C 2 H 4 


C 3 r 6 t PH 3 

C 4 H 8 +H ' 


C 5 H 10 


(13) 


(14) 


in analogy to reactions (7) and ( 8 ) the branching ratio of reaction (13) 
cannot be determined, because- the C 4 H 8 + ion on mass number 56 coincides with 
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the butene ion formed in reaction (9). It appears, however, from the 
pressure dependence of the M56 intensity that the contribution of reaction 
(13) to the mass 56 peak is almost negligible, so that the C3H6'*' forma- 
tion in reaction (13) .must be the predominant one* 

The M29 ion could be conceivably formed from an ethane impurity in 
ethylene, but since the M30 peak was absent, this possibility is excluded. 
The nature of the C2H5"* - ion indicates that it is a product ion, rather than 
a primary ion. However, the pressure dependence of the sum of the inten- 
sities of the M29 sequence is fairly linear as shown in Figure 24, whereas” 
product intensities should show a quadratic pressure dependence around the 
origin of the plot. On the other hand, it may be recalled that signals due 
to this ion sequence, albeit small, were also detected in the 1087A photo- 
ionization of ethylene. Hence, it is reasonable to assume that the C2H 5 + 
ion and the follow-on ion sequence is due to the process 

C 2 H 4 + C 2 H 4 T C 2 H 5 + C 2 H 3 ' (i5) 

The clearer observation of C2H5 + and its reaction sequence at a wavelength 
of 923A compared with 1087$, is mainly due to the higher intensity avail- 
able with the light source combined with a 3 times higher ionization cross 
section of ethylene at 923A. An additional factor may be the participation 
of excited ethylene ions, which according to a mass spectrometric study of 
Botter, et al (Ref. 24) can be formed at wavelengths below about IOI08. At 
923A as much as 2/3 of the total ions produced by photoionization may be 
in an excited state. The quantitative evaluation of the data shows that 
the ethylene parent ion reactions (7) and (15) have the following probabil- 
ities: 83 percent produce C2H5 + ; 11 percent produce G^Hy + ; and 6 percent 
produce C2H5+. The product channel ratio for reaction (7) determined in 
this case is in reasonable agreement with -that found for the same reaction 
in the 1087A photo ionization of ethylene. Any difference beyond the exper- 
imental error must be attributed to the participation of excited ethylene 
ions . 

The formation of butene ions C^Hg*^ via reaction (9) is shown in Figure 
22 in comparison to that determined at 108 7&. As in that case, a linear 
pressure dependence is observed for the intensity ratio at M56 to that of 
all other C2Hg + product ions. However, the rate of M56 formation is only 
about half that observed at 1087^. Again, this behavior is indicative * 
of excited ethylene ions. Due to the additional energy content, the C^Hg** - 
collision complex would have a shorter lifetime, so that its stabilization 
by third body collision is less efficient. 

The origin of the M27 may now be discussed. This ion is not believed 
to be formed from 02^* in a manner similar to M29, because there is no 
evidence for M27 of the follow-on ions 40 and 54 in the photoionization of 
ethylene at 1087A. However, in a mixture of acetylene and ethylene at the 
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same wavelength very weak signals were observed at the corresponding mass 
numbers, indicating the reaction 

G 2 H 2 + C 2 H 4 C 2 H 3 * C 2 H 3 ( 16 > 

Since the ionizing energy at 923$ lies only 0.3 eV above the threshold of 
G 2 H 2 + formation from ethylene, the acetyiene ion cannot carry much internal 
energy and should behave similar to that formed from acetylene at 1087$. 
However, at 923$ the formation of C2H3+ according to reaction (16) would be 
more pronounced due to the higher ion intensity, and because of the absence 
of acetylene which diverts a considerable portion of acetylene ions to the 
M50 and M51 products due to reaction (4). Figure 24 shows also a plot of 
the sum of M27 and follow-on ion intensity versus pressure. The behavior 
is similar to that of the sum of M29 and follow-on ions, i. e. the pressure 
dependence is more linear, than one would expect for a secondary ion sequence. 
It should be noted that in the case of the C2H3 + ion (but not for the 02^"*" 
ion), there may be a small contribution to its intensity from the dissocia- 
tive photoionization of ethylene. The threshold for C 2 H3 + formation by that 
process lies at 898$, i. e. 25$ toward lower wavelengths. While the applied 
spectral resolution excludes the direct admission of this radiation to the 
ion chamber, radiation of shorter wavelengths may nevertheless be present 
as scattered light. 

The C2H3'*' ion sequence constitutes about 2 percent of the total ion 
intensity, and about 20 percent of the acetylene ion intensity if it is 
counted fully as being due to reaction (16). The ratio of acetylene to 
ethylene ions thus produced by photoionization is R = (C2H2 + ) 0 /( G 2 H 6 + ) 0 = 

0.14, if the M27 sequence is added -to that of M26, and R = 0.12 if it is 
not. These ratios are smaller than those found in other studies of the 
photo ionization of ethylene at 923$, but they are still within the same 
order of magnitude. From the data of Botter, et al (Ref. 24) one obtains 
R = 0.21, the data of Brehm (Ref. 25) yield R = 0.16 and those cf Schoen 
(Ref. 26) R - 0.20. A portion of the discrepancy may be due to the varying 
efficiencies of the electron multiplier detectors used in these studies. 

From electron impact ionization studies of ion reactions in ethylene it has 
been concluded (Ref. 23) that charge transfer of 02^**" to ethylene produces 
some losses of acetylene ions. Since the data obtained for acetylene-ethylene 
mixtures at 1087$ gave no evidence for charge transfer, and the probability 
for this process was shown to be less than 15 percent of the total reaction 
probability, the effect of charge transfer to the (C2H2 + ) 0 /( G 2 K 6 + ) 0 tatio 
must be slight. 

If the reactions assigned to the acetylene ions are reactions (10) and 
(16), the branching ratio found for reaction (10) from the photoionization 
of ethylene at 923$ is C3H3*tyC4H5^* - 2.1. This agrees within experimental 
error with the ratio found for the same reaction in a mixture of acetylene 
and 'ethylene, == 2.4, when the photoionizing wavelength was 1087$. 
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Because of the contribution of a dissociative photoionization process to 
the M27 ion sequence when ethylene is photoionized at 1087A, the efficiency 
of reaction (16) can only be estimated. If half of the M27 sequence ion 
intensity is due to reaction (16), its rate is ten percent of the rate 
reaction ( 10 ), i.e. about 14 percent of C 2 H 2 + ions reacting with ethylene 
yield C 2 H 3 '*". 

Having assigned the various products to reactions of either ethylene 
or acetylene ions reacting with ethylene, the associated reaction rate 
coefficients can be determined from the decrease of the M28 and M26 ion 
intensities when compared to the total ion intensities in the corresponding 
ion sequences. The absolute values of the rate constants thus determined 
are k 4 -fl 2 = 0.88 x 10“^ cc/molecule sec, and ^ 7^.13 = 1. 31 x 10 cc/molecule 
sec. The former value. constitutes only 70 percent of the rate coefficients 
for ethylene ions formed at 1087$. This behavior is presumably associated 
with excited ethylene ions formed at 923$. The latter rate coefficient is 
in excellent agreement with the value found for acetylene ions formed at 
1087$, and thus indicates that, in accord with expectation, the C 2 H 2 + ion 
formed from ethylene at 923$ is essentially identical with that formed at 
1087$ by the photoionization of acetylene. 

Further data have been obtained for reactions of acetylene and ethylene 
ions with methane, and for some reactions involving methane' icns. These 
results have not been evaluated and will be discussed in future reports. 

“Yield of Ions Resulting from Photoionization," A variety of small 
molecules has been studied with the photoionization mass spectrometer to 
determine .the yield of - fragment ions resulting from photo ionization. Of 
particular interest in this work was the investigation of the wavelength 
dependence of product ions near their thresholds of appearance, because it 
allows the determination of appearance potentials, which in turn provide 
information on important thermodynamic properties of the ions and their 
parent molecules. 

The techniques involved have been reported previously (Refs. 27,28) 
and only a brief summary is required here. The photoionization mass. spec- 
trometer employed involves a. 1/2 meter Seya grating monochromator with the 
slits set to provide a spectral resolution of about 3.5A at half-width. 

The monochromator exit slit also serves as the entrance slit to the ion 
source chamber. After passing the ion chamber, the radiation is registered 
on a sodium-salicylate-coated photomultiplier. From the ion chamber, ions 
are withdrawn at right angles to the light beam through an aperture of 
3 mm diameter. A concave repeller provides a concentric field configuration 
which enhances the collection of ions. Subsequently, the ions are focused 
onto the entrance aperture of a 180 degree magnetic analyzer, behind the 
exit orifice of which the ions are registered by means of an electron mul- 
tiplier detector. An electrometer readout circuit and recording equipment 
complete the experimental arrangement. 
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The experimental procedure entailed taking simultaneous recordings as 
a function of wavelength of the photomultiplier current and the ion current. 
For this purpose, the magnet current of the mass spectrometer was a juste 
to allow the desired ionic species to pass the mass analyzer, and the wave- 
length drive of the monochromator was actuated* Two types of light source 
were used which gave overlapping spectra in the wavelength region of inter- 
est: (a) a dc cold cathode type source operated with hydrogen as the filler 
gas provided a usefully intense spectrum between 900 and 1300A; (b) a 
Weissler type repetitive spark source operated with argon gave a uselu y 
intense spectrum between 450 and lOOOX. Both light sources emit line spectra 
but the population of emission lines made available is sufficient y ° 

provide data points not more than I0A apart over most of the investigated 
wavelength range. The ion yield curves thus obtained are sufficiently 
well defined to allow their extrapolation to threshold. Since, in addition, 
polyatomic molecules do not usually exhibit narrow structure in absorption 
and photoionization cross sections, the many line radiation source emp oye 
here suffice for this study. 


To date, detailed data have been obtained for methane, methanol, for- 
maldehyde and formic acid. These substances were introduced to the ion 
source chamber via an appropriate leak valve at a pressure of about 10 lorr. 
At these high source pressures satellite mass peaks due to secondary ion 
reactions were observed but their intensity was insignificant. ^ T ese 
secondary peaks, for example M33 from methanol, could be identifie easi y 
from the fact that their thresholds coincided with that of the parent peaks. 
Although the effect of the satellite ion peaks could be reduced by working 
at lower pressures, it was advantageous to make use of the generally goo 
ion intensities provided in the 1 x. 10-3 Torr pressure regime and to toler- 
ate the appearance of secondary ions* 


Figure 25 shows, as an example, the ion yields for and CH 3 observed 

in the photoionization of methane in the 800-100G& wavelength region. From 
the yield curves the onsets can be determined by extrapolation. We find 
the thresholds for C S 4 + at 976& and for CH 3 * at 874A. After correction 
for the involved spectral resolution, the resulting threshold wavelengths 
and associated energies are 


cat 

4 

973 + 3& 

CH+ 

871 + 3& 


12.74 + 0.04 eV 
14'. 23 + 0*06 eV 


These values are in good agreement with those determined by Dibeler, et _a l < 
(Ref. 29) and by Chupka (Ref. 30)* 

- In a similar manner,' thresholds for various ions have been determined 
for methanol and formaldehyde. These are given in Table 1* The relative 
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TABLE 1 


Photoionization Thresholds for the Appearance of Ions in Methanol and in 

Forma Idehyde 


Threshold Wavelength (X) 

-Compound Ion Observed Corrected Energy (eV) 


ch 3 oh 


ch 2 o 


M3 2 

1145 

± 3 

1142 + 3 

10.85 + 0.03 

"431 

1077 

± 3 ' 

1074 + 3 

11.54 + 0.03 

m 

1033 

+ 10 

1030 + 10 ' 

12.04 + 0.12 

m 

952 

+ 7 

949 + 7 

13.06 + 0.10 

as 

900 

+ 3 

897 ± 3 

13.82 ± 0.04 

414 

883 

+ 5 

880 ± 5 

14.09 + 0.08 


M30 

im ± 2 

1139 ± Z 

10.89 

CM 

O 

o 

+1 

M29 

1046 + 3 

1043 ± 3 

11.89 

± °* 03 

M 2 

807 J 3 

804 + 3 

15.42 

± °’ 06 

•M 1 

715 + 3 

712 ± 3 

17.41 

± 0.07 
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ion yields as a function of wavelength for these compounds and for formic 
acid have not yet been evaluated in detail; but as supplementary informa- 
tion we give in Table 2 the ion distribution for all the compounds inves- 
tigated at the photoionizing wavelength of 801&. 

The threshold energies assembled in Table 1 will now be used to derive 
some important selected thermodynamic properties. From the appearance po- 
tentials of the parent ions of methanol and formaldehyde one' obtains the 
ionization energies 10.85 eV from CH3OH and 10.89 eV for CB^O- These values 
are in excellent agreement with those determined by Watanabe (Ref. 31) 
and collaborators using a different photoionization method. Next consider 
the appearance potentials for CH3**" and CH2* from methanol. After subtrac- 
tion of the spectroscopic ionization potentials (Ref. 32) for CH3 and CH2 
(CH3) = 9.843 eV and I (CH2) - 10.396 eV ; respectively, one obtains the 
dissociation energies for the processes 


CH 3 0H .-*■ CH 3 + OH . 

- ch 2 + h 2 o . i> 2 * 

The values thus derived: = 3.98 + 0.04 and D2 = 3.69 + 0.08 can be util- 

ized to derive the heats of formation of the CH3 and CH2 radicals by virtue 
of the relations: 

= AH f <CH 3 ) + ..\H f (OH) - AH f (CH 3 0H) 

D 2 * AH f (CH 2 ) + AH f (H 2 0) - m £ (CH 3 0H) 

The resulting heats of formation are AHf (CHo) =.1.49 + 0.04 eV = 34.3 + 1 
kcal/mole, and AHf (CH2) — 4. 07 + 0.08 = 93.6 + 2 kcal/mole. These values 
are in good agreement with those derived by Chupka (Ref. 30) and by Dibeler, 
et al . (Ref. 29) from the photoionization of methane and CH3 radicals, al- 
though Chupka 1 s values are more precise. From the threshold of CHO'*' appear- 
ance from CH3OH, the heat of formation of CH0 + can be .derived via the relation 

A (CH0 + ) =-* AH f (H) + AH f (CH0 + ) - AH f (CH 3 0H) 


which assumes the the process responsible for the formation of CHO**’ is 


CH 3 0H + hv -* CH0 + + H 2 +' H 
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TABLE 2 


ION YIELDS FOR METHANE, METHANOL, FORMALDEHYDE & FORMIC ACID AT 80lS 


Substance 


Methane 


Methanol 


Forma ldehy< 


Relative 


m/e 

Ion Yield 

16 

0. 585 

15 

0.414 


32 

0.443 

31 

0.492 

30 

0.009 

29 

0.028 

15 

0.031 


30 

0.433 

29 

0.550 

28 

0.016 


46 

0.462 

45 

0.286 

41 

0.003 

32 

0.028 

31 

0.055 

30 

’.0.003 

29 

0. 100 

18 

0.064 
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Since the threshold of CHO* occurs at a higher energy than that of CH 2 0* 
of CH3O* it is reasonable to assume that the above reaction requires two 
steps, with one of these ions as an intermediate. For example, 

CH 3 OH + hv (CH 3 0 + ) + H 

+ * + 

(CH 3 0 ) -> CHO + H 2 

Here, the asterisk indicates that. the intermediate ion must be internally 
excited to dissociate further in the process of CHO* formation. Since the 
dissociation process may require an activation energy, the threshold of 
CHO* formation can provide only an upper limit to the heat of formation of 
CHO+. The present data yield AHf (CHO"*") < 8.62 + 0.10 eV = 198 + 2 kcal. 
This value is in excellent agreement with the heat of CHO* formation pre- 
viously derived by us from the photoionization of formaldehyde (see below), 
forpiic acid, and acetaldehyde (Ref. 33). It thus appears that the activa- 
tion energy ' involved in the dissociation of the intermediate is negligible. 
From the threshold of CH2O* format ion we find that the energy involved is 
by 0.19 eV of 4.4 kcal higher than predicted on account of the process 

CH 3 OH + hv -* CH 2 0 + + H 2 

for which all the required thermodynamic values are known. Thus it appears 
that the intermediate for CHO* formation is CH3O* rather than CH2O*. On the 
other hand, the heat of formation for CH3O* derived from the threshold of 
its observation is AHf (CH3O*) = 7 . 20 + 0. 03 eV = 166 + 1 kcal, in good 
agreement with a value recently derived by Haney and Franklin (Ref. 34) 
from electron impact mass spectrometry with consideration of the excess 
kinetic energy of- the fragments. 

Similar considerations apply to the data for formaldehyde. From the 
threshold of CHO* formation one obtains AHf (CHO*) = 8.51 + 0.03 eV = 196 
+ 1 kcal. This result has been reported previously (Ref. 33). The forma- 
tion of H2* is due to the process 

CH 2 0 + hv -* CO + H* 

The threshold of 15.42 eV is in good agreement with that calculated for 
this process from known available thermochemical values A (H2*) = 15.40 eV. 
The formation of H* is associated with 

CH 2 0 + hv -» HC0 + H + 

This process allows the derivation of the heat of formation of CHO. Sub- 
tracting from the threshold energy for H* formation the ionization poten- . 
tial of the hydrogen atom, I (H) = 13.595 eV, gives the bond dissociation 
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energy for formaldehyde D (CHO - H) = 3.82 eV = 88 + 2 kcal/mole. The 
corresponding value for the heat of formation of CHO is AHf (HCO) = 0.43 
+ 0.07 eV 2J 10 kcal/mole. From this value, the dissociation energy of the 
HCO radical can be calculated as D (H - CO) = 15.7 kcal. In deriving 
these values the new heat of formation of formaldehyde (Ref. 35), ASf (CH2O) 
= -25.9 kcal/mole, has been used. The value for D (CHO - H) derived here is 
the high value of the two possibilities discussed in the literature. Sim- 
ilarly, the value for D (H - CO) derived here is the low value for the 
discussed two possibilities. Independent evidence for D (CHO - H) = 88 
kcal/mole has been derived previously by Haney and Franklin (Ref. 34). From 
a detailed discussion of previous data, Kerr (Ref. 36) has also favored 
this value. 

Although an analysis of the formic acid data is even less complete 
than the ones discussed above, a check has been .made on the heat of forma- 
tion of the CHO radical by the observation of the threshold for the process 

HC00H + hv HCO + OH + 

which was observed to occur at 69oS, corresponding to 17.97 eV. The value 
for the heat of formation of HCO derived from this process is AHf (CHO) = 
0.45 eV s 10.5 kcal. This is in good agreement with the value found above 
from formladehyde, supporting the high value for AHf (CHO) and for D (CHO-H). 

"Measurement of Ionic Mobilities." The effort under this heading repre- 
sents a continuation of previous work concerning the title subject. It has 
been demonstrated previously (Ref. 37) that the photoionization mass spec- 
trometer technique can be employed to obtain reliable drift velocity data 
for mass analyzed ions produced by photoionization. In this manner it is 
possible, in principle, to preselect the energy state of the ion by an 
appropriate choice of varying wavelengths, and to identify the ion subjected 
to measurement by the employed mass spectrometer. 

The experimental set up and the applied experimental procedures have 
been described in detail previously (Refs. 37,38) and will be summarized 
only briefly here. The cylindrical ion chamber configuration employed for 
this application of the instrument is conventional in that ions formed 
along the axis of this light beam are driven to the ion exit aperture by 
an electric field established perpendicular to the direction of the light 
beam. The plane repeller is supported by a hollow stem so that the gas 
pressure in the source can be measured directly with an attached McLeod 
manometer. As a consequence of the ion source geometry the electric field 
is inhomogeneous everywhere except in the vicinity of the axis of ion re- 
moval. The field in this region was calculated from the applied voltage 
and from the distance between the repeller and the ion exit aperture, the 
positions of which were adjusted to a spacing of 0.70 cm. A small correction 
was required to take into account the slight field inhomogeneity near the 
repeller. This was shown by a field plotting method previously (Ref. 37). 
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The light, source is a Weissler type repetitive spark source, providing 
120 pulses per minute of about one microsecond duration each. This time is 
sufficiently short to provide an essentially instantaneous ion deposition, 
unless the residence time of the ions in the ion chamber becomes comparable. 
Under the applied conditions of pressure and electric field, ion residence 
times are usually longer than 10 microseconds. A further necessary experi- 
mental condition for mobility studies is that the mean free path of the 
ions is short compared with the drift distance. This condition is met at 
pressures above about 80 microns. 


Due to diffusion, the ion transient broadens while drifting from the 
plane of origin toward the ion exit plane. Thus, the initially rectangular 
ion transient profile is changed to a Gaussian profile. The average ion 
residence time, t q is defined as the time span required for the center of 
the profiie to reach the exit aperture after deposition of the ions by the 
Accordingly, the ion drift velocity is v - d/T 0 , with d being 
^ e ^^ ft . d:LS 5 anCe < ' 0 * 3 cm in the P resent arrangement). Due to the process 
ot diffusion broadening, some ions arrive at the exit at an earlier time 
t < t Q . It has been shown previously that t and are related bv the 
equation . J 


- 1 - (f/d) (2 Dt) 


1/2 


where D is the ionic diffusion coefficient and f is a factor depending upon 
the threshoid sensitivity for the detection of the onset of the ion pulse, 
f^the detection limit is 10 percent of the total integrated ion pulse, 
f = 1*28, it is five percent of the integrated ion pulse, f = 1.65. 

The former factor has been used previously (Ref. 38), but due to the in- 
creased sensitivity, the latter value will be used here. Then, the measure- 
ment of both t and 1 0 will provide the ionic diffusion .coefficient 


0 = 1 c 
2 f 2 


(1 - t/t Y 


= 1.65 x 10 


-2 


(1 - r/r Y 


Knowledge of the diffusion coefficient and the ionic mobility \i = v/E pro- 
vides the average energy content of the drifting ions. To- derive the random 
energy of the ions, expressed as the ionic temperature, T ± , we make use of 
a relation due to Einstein (see for example Huxley and Crompton, Ref. 39). 


where e is the ionic charge, k is the Boltzmann constant, and F is a factor 
depending on the velocity distribution function of the ions. For a Ifexwellian 
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distribution F - 1, for other distributions F is usually somewhat smaller 
than unity. Here it will be assumed that F » 1. Then the ion temperature 
is defined as 



1 

2 


d 

f 2 


(1 - t/t o ) 
t/? 


■ 6.38 x 10 2 


a - t/t o ) 2 

—nr 


E 


The formula then reduces the ion drift velocity, the ionic diffusion co- 
efficient, and the ion temperature to the residence times t and i Q in the 
ion chamber. These quantities .are determined experimentally from the time 
delay of the ion pulse arriving at the mass spectrometer detector, after a 
group of ions is formed by the light pulse. Delay times are measured as 
a function of the repeller potential and a graphical extrapolation toward 
infinite repeller potential is applied to enable a separation of the ion 
source residence times from the mass spectrometer flight time of the ions. 
The details of this procedure have been published (Ref. 38). 


Below we report data on the parent ions in three gases: argon, nitro- 
gen, and oxygen. 

Argon: The drift velocities observed for A + in argon are plotted 

in Figure 26 as a function of the reduced field strength E/P 0 , for three 
values of the electric field and various settings of pressure in the range 
70-210 microns. The lower group of data points represents measurements with 
the monochromator set to transmit a wavelength of 786& (with 3 . 5 S halfwidth), 
so that only argon ions in the ^ 3/2 ground state are formed. The upper 
group of data points was obtained at 767X. In this case excited %?£/2 argon 
ions are formed in addition to those In the ground state. At 786^, the 
ion transient observed with the oscilloscope was farily symmetric and the 
center of the pulse was used to derive the drift velocities. ■ At 767& an 
additional component was present, the ion transient had an asymmetric appear- 
ance, and the peak of the pulse (which was used to derive the drift vel- 
ocity) shifted forward. The two types of oscilloscope traces are shown in 
Figure 27 for comparison. Clearly, the excited argon ion moves faster than 
the ground state ion. Since the mobility of an ion, like Its diffusion co- 
efficient, • is inversely proportional, to its cross section for momentum 
transfer to the neutral gas; suitably averaged over the range of velocities, 
the data indicate that the momentum transfer cross section for the excited 
argon ion is smaller than that for the ground state argon ion. For both 
types of argon ions moving in the parent gas, the momentum transfer cross 
section is governed mainly by charge transfer. Theoretically one would 
predict, on the assumption that the reduced matrix elements governing the 
two charge transfer cross sections are nearly equivalent, that the differ- 
ence in cross sections is due to the differences in the degeneracy of the 
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Figure 27, 


Argon ion transients* 
upper trace: 76 ?S* 


Lower trace : 786X 
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two processes. Thus* the summation over all possible magnetic quantum num- 
ber combinations indicates Qg^oimd /Q^xcited = because the resonant 
process A+ ( 2 Ri/ 2) + A (^S 0 ) is 2 fold degenerate. The observed increase 
in mobility when going from the ground state to the excited state is about 
three at the lower end of the reduced field strengths. Towards higher 
fields, the ratio becomes smaller due to the fact that the drift velocity 
for the excited argon ions increases with E/P approximately with 1/2 power, 
whereas the drift velocity of ground state argon ions increases essentially 
linearly with E/P* 

The data obtained for ground state argon ions are found to be in 
excellent agreement with previous measurement. Two sets of previous data 
are shown in Figure 26, those by Blond i and Chanin (Ref. 40) are represented 
by the solid line and those by Beaty, e t al . (Ref, 41) are indicated by the 
dashed line, TWo more recent measurements, like the present one with mass 
identification, were made by McAfee, et al. (Ref. 42) and by Mad&on and 
Oskam (Ref. 43). They are in good agreement with the former and the pre- 
sent measurements. An extrapolation of the present data to zero electric 
field leads to a reduced zero field mobility of A + of 4 = 1.5 cm 2 /sec volt. 

It must be realized that the drift velocities of the excited argon ions 
cannot be measured with the same reliance as that of the ground state ions, 
because the two ion transients cannot be sufficiently separated so that 
the arrival time of the center of the excited ion transient is falsified by 
the contribution due to the normal argon ions. Accordingly, the data pre- 
sented here for that ion species must be considered tentative* Qualitatively, 
however, there can be no doubt about the increased drift velocity of the 
excited ion compared with the ground state ion. The importance of this 
finding should not be underrated. It has been made possible solely by the 
energy differentiation provided by the method of photoionization. 

The good results obtained for ground state ions in argon have en- 
couraged measurement of diffusion coefficients and ion temperatures accord- 
ing to the methods outlined above. These are shown in Figures 28 and 29. 

The straight line in Figure 10 is calculated from the thermal, zero field 
mobility of the argon ions. The data points agree with the calculated 
diffusion coefficients, indicating that the energy that the argon ions 
acquire from the field is quickly passed on to the neutral argon matrix. 

This behavior is reflected also in the derived ion temperature, which lias 
only little above that corresponding to thermal ions* A considerable 
scatter of the data points in Figure 29 makes apparent that an accurate 
determination of ion temperatures Is not possible by this method* However, 
the good agreement of ion temperature and diffusion coefficients with the 
expected ion behavior indicates the validity of the concepts underlying 
their derivation. 

Nitrogen: Drift velocities observed for nitrogen ions formed from mol- 
ecular nitrogen at an ionizing wavelength of 790A were reported previously 
by us and are shown in Figure 30* There is good agreement with the results 
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Figure 28; Diffusion coefficients for A + ( 2 P ) ions in argon, 

line calculated for thermal ions. 
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Figure 30. Drift velocities of nJ in nitrogen versus relative field 
oJ^points 7648 iZing Wavelength: lower grou P of Points 79oS; upper group 



of Saparoschenko (Ref. 44) and Keller , et al . (Ref. 45) who also used mass 
identification of the ions. Similar to the behavior of argon ions, an 
increase of drift velocities was observed for N2 + as the ionizing wave- 
length was lowered. However, the mobility increase occurred at isolated 
wavelengths, and was then clearly recognizable as an asymmetry of the ion 
transient. A very strong asymmetry occurred at 764$, and drift velocities 
corresponding to the peak of the ion pulse at this wavelength are shown in 
Figure 30 as the upper group of data points. The increase in ion mobility 
appears to occur at wavelengths where the emission lines available from 
the light soure coincide with autoionization features in the N£ absorption 
spectrum. At 764$, the wavelength investigated here, the line group emitted 
from the source overlaps the v (1-0) transition at 764. 35$ of the Ogawa 
and Tanaka (Ref. 46) progression (1) series which, having a high ionization 
efficiency produces a considerable population of vibrations lly exicted 
N2’*" (v = 1) ions. Accordingly, the increased mobility, v, at 764$ is 
due to vibrationally exicted N2"** i° n s, whereas at 790$ only ground state 
nitrogen ions are formed. Similar to the results discussed above for 
argon ions, the N2* ion mobility increases by a factor of about 3 for 
vibrationally excited ions. A similar behavior with E/P is observed as 
for ions in argon. 

Figures 31 and 32 show the diffusion coefficients and the ion tem- 
peratures determined for N2* ( v ■ 0), respectively. Unlike' the behavior of 
argon ions, the diffusion coefficients are much higher than that for ther- 
mal ions (indicated by the solid line), and the ion temperatures derived 
lie also considerably above thermal values. Both factors indicate an 
appreciable energy content of N2* ions in the presence of electric fields. 
Note also that both .diffusion coefficients and. ion temperatures increase 
as the field strength is increased. 

Oxygen: The drift velocities determined for 02 + ions as a function of 

E/P are shown in Figure 33. These data were obtained with an ionizing 
wavelength of 1007$, and good agreement is observed with the results of 
previous investigations (Refs. 47,48) in which the identity of the ion was 
inferred rather than determined. At 1007$, photo ionization of O2 may 
lead to 02* in the first vibrationally excited state, but from the results 
it appears that the direct photoionization probability for the formation 
of 02^ (v = 1) is only a very small fraction of that leading to ground state 
02“*“ ions. Unfortunately, a sufficiently strong emission line was not avail- 
able in the region 1007-1027$ (where O2 ionization initiates), so that the 
O2"*' (v = 1) formation could not be positively excluded. However, similar 
to the previous case (N2 - **), an increase in mobility was observed at auto- 
ionization features of oxygen. Figure 34 shows -that wavelength at which 
the 02 + residence time in the ion source chamber decreased when the other 
experimental conditions are kept constant. Each such decrease of residence 
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Figure 31: Diffusion coefficients for N 2 (v -= 0) ions in nitrogen. 

Solid line calculated for thermal ions. 
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Figure 34: 
tion features 
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Correlation of oxygen on residence time with autoioniza- 
in the absorption spectrum of oxygen. 
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time (or increase of ion mobility) correlates with an autoionization feature 
in the absorption spectrum of oxygen in the 900-1000& wavelength region 
here investigated. The limited number of strong emission lines available 
from the light source prevents the demonstration that an increase of mo- 
bility is associated with every sufficiently energetic autoionizing feature, 
but the number of available lines appears to be sufficient to prove the 
point. The variation in residence times observed is due to the degree of 
overlap of emission and absorption intensities and the varying ionization 
efficiencies. 

Figures 35 and 36 show the diffusion coefficients and ion temperatures 
derived for C>2 + formed at 1007 a. Again it is evident from both groups of 
data that the ions drifting in the electric field are non -thermal, although 
the increases in the diffusion coefficient and the ion temperature above 
thermal values are not as drastic as was found for nitrogen ions. 

B. Theoretical Studies 


Published Data . - 

"Mariner 6: Origin of Mars Ionized Carbon Dioxide Ultraviolet 

Spectrum" 

A. Dalgarno, T. Degges and A. Stewart 
Science 167, #3924, 1490 (1970) 

Predicted intensities of the ionized carbon dioxide (C02 + ) emission 
feature at 2890A and the Fox-Duffendack-Barker bands are 5.2 and 19.9 kilo- 
ray leighs, respectively, for a vertical column. Direct photo ionization 
of carbon dioxide by solar radiation contributes 3.5 and 4.1 kilorayleighs, 
respectively, and fluorescent scattering by C0 2 + , 1.6 and 15.3 kilo- 
rayleighs, respectively. Photoelectron impacts are less important. 

,r Dipole Properties of Molecular Hydrogen" 

G. Victor & A. Dalgarno 

Journal of Chemical Physics 50, #6, 2535 (1969) 

Experimental data on electron scattering, optical refract ivity, and - 
molecular anisotropy are used to construct a model dipole spectrum for 
molecular hydrogen which is consistent with oscillator strength sum rules. 

The model spectrum is used to calculate the parallel and perpendicular dynamic 
dipole polarizibilities, the refractive index, the Verdet constant, the 
Rayleigh scattering cross section, the Rayleigh depolarization factor, and 
the molecular anisotropy as a function of frequency of the incident rad- 
iation. The tensor components of the van der Waals coefficient for a pair 
of interacting hydrogen molecules are also calculated. 
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Figure 35: 


Diffusion coefficients for 0^" ions in oxygen. Solid line 
— i- i for thermal ions. 
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"The Effect of Oxygen Cooling in Ionospheric Electron Temperatures" 

A. Dalgarno, M. McElroy, M. Rees, & J. C. G. Walker 

There are discrepancies between ionospheric electron temperatures de- 
rived from Thomson scatter data and electron temperatures predicted from the 
solar ultra-violet heat source. The inclusion of electron cooling by 
excitation of the fine structure levels of atomic oxygen removes the dis- 
crepancy ' throughout the day at all altitudes above 320 km. Belov; this al- 
titude a discrepancy persists, but it probably lies within the uncertainties 
arising from the basic atomic and molecular data employed in the theoretical 
analysis and from the solar flux data. 

"CO* Day glow on Mars and Venus" 

A. Dalgarno and T. Degges 

Quantitative calculations are presented of the intensities of the 
C02 + bands present in the dayglows of Mars and Venus and it is argued that 
fluorescent scattering by CO 2 * icns and photoionization of CO 2 are the main 
sources of excitation. An estimate is made of the intensity of 3914$ 
emission that would arise if the atmosphere contained N2* 

"Electron Impact Excitation of the Dayglow" 

A. Dalgarno, M. McElroy and A. Stewart 
Journal of Atmospheric Sciences 26, #4, 753(1969) 

Calculations are described of the equilibrium velocity distributions 
of the photoelectrons produced in the F -region by solar ionizing radiation. 
Detailed estimates are presented of the intensities and altitude profiles 
of emission features of atomic oxygen, molecular nitrogen, and molecular 
oxygen, appearing in the dayglow as a result of photoelectron impacts. 

"Carbon Atoms in the Upper Atmosphere of Venus" 

F. F. Marmo and A. Engelman 

ICARUS 

The detection of carbon atoms in the upper atmosphere of Venus is re- 
ported on the basis of an extension of a previously reported analysis of 
low resolution spectra of the solar illuminated atmosphere of Venus. It is 
demonstrated that the previously reported anomalously high albedo for W > 
1500$ can be explained by invoking the presence of a 14 kR Cl signal at 
1657$. It is possible to check the validity of this conclusion on a future 
Venus f,ly-by mission which employs a scanning spectrometer operating in 
this spectral region. Furthermore, if it is valid to extrapolate these 
results to the case of Mars the vacuum ultraviolet spectra acquired by 
the Mariner 6 and 7 -missions should include this spectral feature as well as 
some of the more intense bands of the CO 4th positive system. 
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Unpublished Data .- ,f Mars Lander Experiment — Spectral Photometric 
Observations of the Martian Atmospheric Dayglow. 11 The identification of 
the atomic and molecular constituents in a planetary atmosphere is a funda- 
mental branch of planetary -aeronomy. For the earth atmosphere., one of 
the most prolific sources for these type of data have been obtained through 
the performance of airglow observations employing ground based spectro- 
photometric techniques. However., the performance of such experiments are 
severely constrained owing to two important limitations: (a) the spectral 
region of investigation is limited to wavelengths greater than 3000$ owing 
to , atmospheric attenuation, and (b) throughout the visible regions intense 
background radiations .(on the order of megaray leighs) can be encountered 
owing to Rayleigh and Mie scattering of the solar illuminated atmosphere. 

The identification of organic and/or inorganic constituents in the Martian 
atmosphere is fundamental to the associated aeronomic and exobio logical 
studies. For this purpose the established techniques utilized from the sur- 
face of the earth can be translated, extended and modified into an appro- 
priate configuration so that airglow Martian surface observations may be 
performed taking full advantage of the reduced background scattering as 
well as the relative absence of atmospheric absorption down to 2000$. Ow- 
ing to the great potential involved a number of relevant initial tasks 
have been performed in order to evaluate the feasibility of performing 
such a program as well as to identify the type and variety of optical data 
that can be acquired from a surface based Martian airglow station. In 
order to operate under realistic constraints it is proposed to employ a 
rugged, double pass, scanning spectrophotometer probing the spectral region 
\\ 2000-8000$. The configuration represents a simplified version of a 
fully developed instrumentation currently used in earth satellite studies. 
The proposed spectrometer involves no moving optical components so that 
it can be packaged for extreme ruggedness in order to survive the rela- 
tively high g-loading in the landing environment. 

The specific experiments suggested in the following sections involve 
a host of complex theoretical and practical problems which are considered 
to lie outside the scope of the present survey and feasibility study. In 
fact, in some cases, it is difficult to assess the value of performing a 
specific experiment as compared to other techniques which may be available 
or have been suggested for the same purpose. Alternatively^ it Is felt 
that enough scope and data have been incorporated into the present dis- 
cussions to amply illustrate the potential involved in a number of applica- 
tions of airglow observation technology for the identification of specific • 
atomic and molecular constituents in the atmosphere of Mars- 

The remaining discussion is presented in the following four sections. 
The pertinent geometric and solar Illumination factors involved in the per- 
formance of the proposed Martian twilight, day and. night airglow experiments 
are discussed in addition to a number of possible signal generation sources 
within the spectral region XX 2000-8000$ involving some of the more possible 
constituents in the Martian atmosphere, consideration of the roles of re- 
sonance scattering, fluorescence scattering, atmospheric absorption. 
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nightglow, and Martian dust, and the experimental configuration of the 
suggested scanning spectrophotometer is discussed. 

Most of the pertinent geometric parameters involved in this study 
are illustrated in Figure 37. In this section, a number of relationships 
will be deduced between these several parameters which are employed sub- 
sequently in' the development and discussions of the' suggested experimental 
configurations . 

Zenith (p=0) sky radiances have been calculated for both earth and 
Mars for the case of solar zenith angle of zero degrees (i.e., overhead 
sun) employing available model atmosphers for earth and Mars, (Refs. 49,50) 
the appropriate incident solar fluxes (Ref. 51) and the absorption (Ref. 52) 
and Rayleigh scattering (Ref. 53) characteristics of the several pertinent 
planetary atmospheric constituents. The results are shown in Figure 38 
wherein the heavy dashed and solid curves obtain for the earth and Mars 
zenith radiances respectively. It should be noted that two separate rad- 
iance scales are employed in this figure which illustrates two important 
factors. First, over the spectral range \\ 3000-8000$, Martian atmospheric 
background intensities are about too orders of magnitude less than those 
on earth (which is roughly equivalent to a platform located at about 30 km 
above the surface of the earth) and second, the fact that the observational 
spectral range is extended down to 2000$ for the case of Mars. However, 
a portion of the apparent improvement in . background radiation is due to 
the scaling factor applied to the lower incident Martian solar flux. Thus, 
solar flux dependent signal sources would undergo a corresponding decrease 
whereas in the case of self-emission sources, full advantage can be taken 
of the lower background radiance conditions. It should be stressed that 
only the role of Rayleigh scattering has been considered in these calcula- 
tions so that the presence of a solar -illuminated surface haze layer would 
result in a Mie scattering background intensity contribution which has not 
been included herein. In any event, it appears that as far as background 
considerations are concerned, distinct advantage is obtained with respect 
to the terrestrial atmospheric case for the performance of dayglow obser- 
vations from the Martian surface. 

For the earth case, it has been well established that twilight air- 
glow observations contribute significant data to the aeronomy of the upper 
atmosphere. In the following discussions, reference is made to the per- 
formance of cogent twilight observations in both the zenith (p = 0) and 
horizon (or more appropriately at j3 RJ 80°) directions. Accordingly, for 
subsequent reference it Is convenient to derive a number of relationships 
involving the geometric parameters identified in Figure 37. 

The twilight geometry pertinent to planetary atmospheres has been dis- 
cussed in detail elsewhere (Ref. 54). For the simplified case where at- 
mospheric refraction is neglected, a sharp diadow line AB prevails as shown 
in Figure 37. Under such conditions, it can be shown that for the zenith 
look-angle, p = 0° and small a-values. 
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for earth 
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z 

2h £ a 2 for Mars (l 8 ) 

z 

In addition, for horizon observations performed along (3 - 90 , the horizon 
illumination altitude, h H as shown in Figure 37 was determined directly 
from geometric relationships. The zenith and horizon illumination altitudes, 
h z and h^ respectively for Mars are shown in Figure 39 in terms of both the 
solar depression angle, a, in degrees, and in minutes. 

Resultant zenith radiances are presented in Figure 40 in terms of 
Rayleighs /$ as a function of wavelength for XX 2000-8000A and a number of 
h 2 -values ranging from 0-100 km. For convenience, the corresponding a-value s 
are included on each h z - curve. The horizon radiance in Rayleighs /A is 
presented in Figure 41 as a function of solar depression angle. In this 
figure, the radiance calculations were performed for three representative 
wavelengths, namely 2000, 2500, and 3000$ since no significant attentuation 
occurs for XX > 3000$ where the screening height approaches zero km. The 
screening height, h s results from attentuation along the twilight path 
indicated in Figure 37 as A’B\ In the case of Mars, significant attenua- 
tion occurs for XX < 1975$ where C0 2 absorption predominates. At higher- 
wavelengths, the decreased attenuation is due only to atmospheric Rayleigh 
scattering by the major constituent model atmosphere gases, i.e. , C0 2 and 
N2« (Ref . 50). .The sharp increase in h s -values for XX < 2000$ shown in 
Figure 42 illustrates this point. In the present discussion, the h s -values 
are small compared to the radius of the planet at all wavelengths XX 2000- 
8000$. Under this condition and for relatively moderate solar depression 
angles, it is appropriate to simply add the h s and h z ~values to derive the 
effective altitude above which resident atmospheric constituents are fully 
solar illuminated under given twilight conditions. This is illustrated 
in Figure 37 by the altitude region above the point marked B 1 . 

Estimates were made for the time duration of specific twilight condi- 
tions which obtain on both earth and Mars for comparison and reference. 

For the present purpose, dn arbitrary twilight condition is defined for 
a change in solar depression angle, £C£ from 1 to 18 degrees (which corres- 
ponds to astronomical twilight on earth). The time durations of this con- 
dition as a function of date at common latitudes were acquired for both 
planets. The data, pertinent to the case of earth was deduced directly 
from data contained in the American Ephemeris and Nautical Almanac. Al- 
though corresponding data for Mars were 'not readily available, they were 
calculated by application of the following expression 


A . | o T , 

At T = 1 360/24.7 » 

At T « time duration for astronomical twilight = Aa - 0-18 , cos t q « 
-tan 5 tan A. 


(.19) 
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Figure41 • Horizonsky radiance as a function of solar depression 
angle for X = 2000, 2500, and 3000A. 

7R 
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Figure 42, Screening height as a function of wavelength for Martian 
twilight. 79 



sin g> sin X 4- sin a 
cos 6 cos X 


( 20 ) 


where 


and 


5 = solar declination angle 
X = geographic latitude 
a ~ solar depression angle 


The factor 360/24.7 relates to the rotation period pertinent to Mars. 

The time duration- of astronomical twilight (Ax = 0-18°) for the northern 
hemisphere of earth and Mars, are shown in Figures 43 and 44 respectively. 

The similarity of the At<£ -values for common latitude is striking although 
this behavior reflects the ensemble of pertinent planetary and orbital 
characteristics of the two planets. 


In summary then, the following pertinent features can be noted: (a) 
for the spectral region XX 3000-8000A, ' the zenith sky radiance observed on 
the surface of Mars is about two orders of magnitude less than that which 
prevails on the earths surface and that this background corresponds to 
that encountered from a platform located at a 30 km earth altitude, (b) 
for Mars, the spectral region of investigation can be extended down to 
2000?, (c) the time duration of specified twilight conditions for earth 
and Mars are similar, and (d) for given a- values, the shadow height on Mars 
is about half of the earth value, which allows increased observation time 
per unit altitude interval. 

The relationships and results presented in Figures 37 through 44 in- 
clusive will be referred to and employed as required in the following dis- 
cussions which identify a number of possible signal sources which emit 
throughout the spectral region XX 2000-8000?. 


Herein a number of possible signal sources are discussed and evaluated 
in order to illustrate the potential involved in performing the proposed 
Martian Lander spectrophotometric experiment. Owing both to the broad 
scope of the problems encountered and to the required brevity of the .pre- 
sentation, detailed discussions involving parameters which are of necessity 
somewhat speculative have been avoided. Thus, each signal source category 
discussed results in the illustration of the potential involved rather than 
in the recommendation of specific detailed selective experiments. A con- 
venient format has been adopted which incorporates five categories: (1) 
Resonance Emission from Solar Illuminated Atoms or Ions in the Upper Atmos- 
phere, (2) Fluorescence Emission from Solar Illuminated Molecules or Mol- 
ecular Ions, (3) Airglow Radiations from Chemiluminescent Reactions between 
Atmospheric Constituents, (4) Absorption Owing to the Presence of Trace 
'Constituents in the Lower Atmosphere of the Planet, and finally, (5) The 
Role of Particulate Scattering due to Dust and Haze in the Lower Atmosphere 
of Mars. 
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In general, atomic, ionic and meta stable species in a planetary atmos- 
phere reside in the upper altitudes owing to recombination processes which 
become more efficient with increasing' number densities. The resultant 
high altitude profiles are particularly amenable to observations conducted 
under twilight conditions wherein the background radiation can be minimized 
without significant reduction in the signal intensity. This 'is not to say 
that full dayglow observations are precluded; thus for high signal intensity 
magnitudes, a program of dayglow observations may be warranted. 

With certain simplifying assumptions. for the' case of resonance scatter- 
ing by atoms or ion species 


, where 4fll 

X 

f 

b) 

N 


2 

4kI - «F ~ X 2 fw - ( 21) 

me ^ 

= the emission rate in photons-cnT^-sec - ^-?"* 1 

6 —2 —1 
(one Rayleigh - 10 photons -era** -see” ) 

- the solar flux in photons-cm”^-sec"^-S**^ 

= resonance wavelength 

* f-value or oscillator strength of the line 
= the albedo for single scattering 

-2 

number of atoms along the line of sight in atoms-cm 
= cos 6 where 0 is the angle between the illuminating and 
observing paths. 


For the present purposes, it is convenient to define an emission probability ' 
factor, I.e, the number of photons scattered per second per solar illumin- 
ated atom: 


so that 


g 



fu 


( 22 ; 


4*1 = g 


(23) 
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The detailed composition of the Martian atmosphere is unknown at pre- 
sent. However, even with the limited available data, one may speculate 
that an ensemble of possibilities exist for the production of neutral and 
metastable atoms and ions in the upper atmosphere. When solar illuminated, 
these species will re-emit at their characteristic signature resonance 
radiations which can represent possible signal sources of interest for W 
2000-8000A. Among the generation processes are included: ionospheric re- 
combination, solar photodissociation, photoelectron and solar wind proton 
interaction with the atmospheric constituents, photochemical production of 
these species, etc. Under certain conditions, it can be expected that the 
resultant resonance signals could be highly intense. For example, in the 
terrestrial atmosphere case the [Ol]2i“6300i£ and the sodium-D-589o£ radia- 
tions (on the order of 20 kR) have been observed in full dayglow experi- 
ments by Noxon and Goody (Ref. 55) and Donahue (Ref. 56) respectively, 
in spite of the intense earth background radiations shown in Figure 38 at 
these spectral positions. Since the corresponding background intensities 
pertinent to Mars are two orders of magnitude lower, it is evident that in 
the proposed Martian surface experiment similar and lesser emission inten- 
sities could be observed in the dayglow. Another possible source of atmos- 
pheric atoms and ions results from the deposition and subsequent observa- 
tion of interplanetary debris in the upper atmosphere of Mars. This latter 
source is singled out for more detailed discussion in order to illustrate 
the potential involved. 

Interplanetary debris represents an important source of minor consti- 
tuent atoms and ions in the upper atmosphere of the earth. The identity 
and inventory of these debris atoms and ions have been experimentally con- 
formed for Na*, Mg***, A#*, Ca+, Fe* and Si* (Ref. 57) by direct in-si tu 
measurements' and also for Ca*, Na, K and Li (Ref. 58) by ground-based 
photometric twilight observations. The observed atomic number density-al- 
titude profiles typically peak in the vicinity of 100 km with a sharp 
cutoff occurring at lower altitudes due to the occurrence of oxidation 
processes. It may be speculated that a somewhat similar situation obtains 
in the upper atmosphere of Mars. In fact, on the basis of an analysis 
performed by Marmo and Brown (Ref. 59), it has been shown that assuming 
equivalent debris influx rates, for the case of Mars the peak number den- 
sity would be larger (due to the limited oxidation processes compared to 
earth) and also it would occur' at altitudes in the vicinity of 80 km. In 
any case, it appears reasonable to search for the presence of these debris 
species in the upper atmosphere of Mars. To emphasize this possibility, 
the following preliminary analysis was performed. 

First, a number of possible debris constituents were selected on 
the basis of too criteria: (a) that they possess a relatively high chem- 
ical abundance in typical meteors (Ref. 60), and (b) that their resonance 
wavelengths' are located at W> 2000& (Ref. 61). The selected debris 
atoms and their corresponding resonance wavelengths are presented in col- 
umns 1 and 2, respectively of Table 3. From readily available data 
(Refs. 61,62) corresponding g-values were calculated as shown in column 3 
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TABLE 3 


. Relative Emission Probability - Abundance Factors (gC-values), for Solar 
Illuminated Meteoric Debris Atomic and/or Ionic Species 
in the Martian Atmosphere 


Debris 

Species 

Resonance 
Wavelength A 

/ photon \* 
s \ sec-atom/ 

gC-values 

A4 I 

3082 

7.6(-2) 

2.6(-l) 

Ca I 

4227 

1.0(-2) 

1.7 (-2) 

Ca II 

3933 

1.4(-1) 


Fe I 

3441 

7. 6 (-3) 

6. 0(-l) 

Fe II 

2599 

2. 2(-2) 


K I 

7699 

1*2(-1) 

1.4(-1) 

• Li I 

6708 

4.4(0) 

1.6 (-2) 

Mg I 

2852 

1.0(-1) 

3.2(0) 

Mg II 

2796 

3.1(-2) 


Na I ‘ 

5890 

6.4(-l) 

1.0(0) 

Ni I 

3370 

8.2(-2) 

7.9(-2) 

Si I 

2514 

6.0(-3) 

2.1(-1) 

Ti I 

3342 

4. 0(-2) 

3. 4 (-2) 


Values in parenthesis are raised to that power of ten 
gC-values are relative to a value of unity for sodium 
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of Table 3. Relative signal intensity estimates can be derived for the 
tabulated species as follows. Firsts it is assumed that the deposited 
inventory of the upper atmospheric debris species is equivalent to its 
chemical abundances in typical meteors. Then if C represents this abun- 
dance parameter , the product gC-values in column 4 of Table 3 are measures 
of the expected relative signal intensities. It is convenient to normalize 
these values to sodium which is assigned an arbitrary value of unity; 

As noted previously, a representative earth value for the Na-D-5890A 
dayglow can be taken as about 20 kR (Ref. 54). If, this signal intensity 
value is employed as a reasonable estimate for the Mars dayglow, it 
follows that all the species indicated in Table 3 should be readily ob- 
servable under certain specified conditions (of both dayglow and twilight 
airglow) encountered during the conduct of the suggested experimental pro- 
gram. 


Molecules and molecular ions which may represent sources of Martian 
atmospheric emission included the following candidates: CO, NO, 02.? N2, OH, 
CN, 0 2 + j C0+, C02 + , and others. Therefore, the resultant fluores- 

cence over the spectral region of interest W 2000-8000^ owing to the in- 
itial absorption of solar radiation by these molecules constitutes an- 
other signal source which is evaluated in the following discussion. 


The molecular case is more complex than its atomic analog discussed 
previously although the single scattering theory for atoms may be adapted 
readily for application to diatomic molecules. The expression which des- 
cribes single molecular photon scattering is analogous to Equation (21) 
for atomic species and is given by: 


4 *W * "Vo *4 £ 

me 


(24) 


where the v'v" subscripts emphasize the fact that the initial absorption 
may occur in one electronic-vibrational transition while the subsequent 
fluorescence may occur at longer wavelengths involving other electronic- 
vibrational transitions. In the above expression, the F, f, and A- values 
have subscripts v’o to indicate involvement of the absorption band of the 
lowest vibrational level (v" = 0) of the ground electronic state. The 
N Q -value indicates that the participating molecules generally reside at 
the lowest vibrational level of the ground state. As a further analogy 
to the case of atoms, the emission probability factor for molecular 
fluorescence may be expressed as 


so that 


te 2 2 

g v'v" = "Vo 2 Vo Vo “v'v" 

me 


4rtl , „ * g , „ 

V*v v V p 


( 23 


( 26 ) 



In general, the pertinent molecular f -values are not available so that 
the corresponding g v * v "-f actors are not readily calculable. However, the 
relative intensity distribution of the molecular resonance and fluorescent 
bands can be deduced for an optically thick atmosphere on the basis of the 
published theoretical results of Chamberlain and Sobouti (Ref. 63). They 
have shown that for this case 


p. 

4*1 , ,, = *F t to , tl AX — H(p ) (& '. 

v’v" v f o v B v" o' 

where H(p. 0 ) are tabulated functions by Chandrasekhar (Ref. 64) and Ax is 
the effective line-width (Ref. 63). Thus, the relative intensities of 
both the resonance and fluorescence bands are proportional to the solar 
flux, F v » 0 and the albedo for single scattering, Uy' v t'« The f v i Q -values. 
are not involved in Equation (27) since the solar flux penetrates into the 
optically thick atmosphere to saturate the bands. If the AX (^ 0 /jji)H- value 
is assigned the value of unity, then 


4*1 * it ^ *F • to a || 
v*v v*o v V 


(28) 


which yields reasonable signal estimates for the optically thick case, as 
well as indicating the relative intensities in a given fully solar illum- 
inated band system. 

A number of pertinent parameters for several possible molecular con- 
stituents in the Martian atmosphere are summarized in Table 4. The band 
systems are identified and the appropriate electronic transitions are 
indicated in column 2 while column 3 includes some selected v'v” transitions 
and their corresponding X-values as presented by Barth (Ref. 65). In 
each case, the initial absorption (i.e. , v 1 = 0) transition can be iden- 
tified easily while the subsequent fluorescent emissions were selected on 
the basis of their relative efficiency over the spectral region XX > 2000A. 
In accordance with Equation (28), the 4*F v t 0 to v .i v >«- values corresponding to 
the selected fluorescent emission bands are presented in column 4 of the 
table. At best, these values represent maximum signal estimates which 
can be expected to be reduced significantly in practice for specific 
cases. Finally, some general comments are presented in column 6 which are 
pertinent to each of the individual transitions. 

Detailed evaluation of the probability of observing any specific 
emission indicated in the table would require a complex analysis where 
much of the basic input data which at this point are at best speculative 
or just not available. For example, it is important to ascertain whether 
a particular molecular constituent is well-mixed or is resident in layer 
form in the upper Martian atmosphere. It can be assumed that the major 
model atmospheric constituents, N 2 and CO 2 are mixed, while for such 
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TABLE 4- 




Absorption- 

F . to - value, 
V ° 2 


Molecular 

Band System 

Fluorescence . 


Constituent 

Transition 

Transition 
v'-v M X 

photons /cm /sec/A 

Co laments 


CO 

4^ Pos. System 
A ^ - X 1 2 + 

5-0 

5-11 

5-12 

5-13 

1393 

2012 

2089 

2172, 

2.8 (7) 
1*5 (7) 

1.8 (7) 
1.1 (7) 

Allowed Transition 

Screening height for 1393$ & 160km 

(see Fig. ) 


Hopfield- 

0-0 

1804 

2.9 (6) 

Forbidden Transition 


Birge Bands 
a' 3 £ + - X‘ X E + 

0-5 

0-6 

0-7 

2223 

2327 

2440 

1*1 (9) 
1*3 (9) 
1.3 (9) 

Screening height for 1804$ ~ 90km 


Cameron Bands 

0-0 

2063 

2.3 (10) 

Forbidden Transition 


a Vi - X l Z* 

0-1 

2158 

2.8 (10) 

Screening height for 2063? % 20km 
Very rich fluorescence spectrum for 
XX 2000-2800X (a) 

NO 

Gamma Bands 

0-0 

2262 

5.9 (10) 

Allowed Transition 


2 e + - x hi 

0-1 

2363 

8.2 (10) 

Screening height for 226 2& X 15km 


1-0 

2149 

9.8 (10) 

Screening height for 2149? # 15km 

° 2 

S chumann-Rung e 

0-0 

2026 

8.0 (2) 

Allowed Transition 


3 - 3 - 

0-12 

3105 

8.0 (9) 

Screening height for 2026X £ 20km 


Z - X 
« 8 

0-13 

3234 

8.8 (9) 

Fluorescence spectrum rich in bands 
for XX 1700X - 4000X .(b) 


Her zb erg Bands 

0-0 

2885 

2.4 (7) 

Forbidden Transition 


2 3 - 

0-8 

4370 

3.8 (11) 

Screening height for 2885A ~ 5km 


A Z -X 

0-9 

4647 

3.9 (11) 

Fluorescence spectrum rich in bands 


u g 

0-10 

4955 

3.2 (11) 

for XX 2400-5000? (c) . 


Molecular 

Constituent 


N 


2 



Absorption- 

F * w - value 

Band System 

Fluorescence' * 

v*o 2 

photons /cm / sec /A 

Transition 

Transition 


vVL 

X 


Vegard -Kaplan 

0-0 

2010 

7.8 (7) 

Bands 

0-4 

2463 

8.5 (9) 


0-5 

2605 

l.Cr (10) 

A 2 £ + - X l T + 

0-6 

2762 

9.7 (9) 

• u g 





Hopfield- 
Birge Barids 


hi 


V 

8 


0-0 


986 


Comments 


Forbidden Transition 

Screening height for 2010A « 20kra 


Allowed transition 

Screening height for 986A > 200km 

No significant fluorescence for 



Lyman-Birge- 
Hopfield Bands 

a H - X l L 
g 8 

0-0 

0-9 

0-10 

0-11 

1450 

2040 

2130 

2227 

1.5 (7) 

5.0 (4) 

7.1 (3) 
7.3 (2) 

Forbidden transition 
Screening height for 1450A a* 150km 
Photon emission flux for X > 2000A 
too weak to observe 

OH 

Hydroxyl UV 
Bands 

A 2 £ + - X hi 

0-0 

1-0 

3064 

2811 

4.0 (12) 
5.4 (11) 

Allowed transition 

The g-values for the transitions 0-0, 
. and 1-0 are 8.6 (-4) and 2.4 (-5) 
photen-sec-atom respectively (d) 

ON 

Cyanogen Violet 
Bands 

B V - X 2 Z + 

0-0 

1-0 

3876 

3582 

8.1 (12) 
9.2 (11) 

Allowed transition 

The g-values for the transitions 0-0, 
and 1-0. are 7.4 (-2) and 7.3 (-4) 
photon-sec-atom respectively (e) 

4 

First Negative 
Bands 

0-0 

3914 

6.4 (12) 

Allowed transition 

From published £ -values, 4 x 10 

ion-cm 2 results in signal intensity 


B 2 y + _ x 2 Z + of 1 kR. 

“u. S 



Molecu lar 

Band System 

Absorption- 

Fluorescence 

F , <o - value 
v'o 2 


Constituent 

Transition 

Transition 

photons/cm /sec/a 

Comments 



vV 

X 



°2 

Second Negative 

0-0 

2610 

8.0 (6) 

Allowed transition 


System 

11-0 

2153 

2.0 (10) 

Screening height for 2610 & 15km 


A hi - X hi 
U - 8 

11-1 

2242 

6.0 (10) 


co* 

First Negative 

0-0 

2191 

1*4 (11) 

Allowed transition 


Band 




Screening height for 2191a # 15km 


2 + 2 + 
b - x ir 




The g-value for the transition 0-0 
is 1. 2 x io** 4 photons-sec-atom (f ) 


Comet Tail Bands 

0-0 

4900 

2.2 (12) 

Allowed transition 


a *n - x 2 z + 

0-1 

5487 

5.5 (12) 

The g-value for the transition 2-0 
is 1.3 x 10" ^ photon-sec-atom (g) 

cot 

Fox, Duffendale 


3378 

5.0 (11)* 

Allowed transition 

4 

and Barker Comet 


3505 

5.0 (11) 

Intense bands seen in comets by 


Bands 


3674 

5.0 (11) 

fluorescence radiation, (h) 


9 9 




Believed to be initially excited by 


2 n - x ji 

U g 




X « 2800&; screening height £ 5km. 


(a) Barth, C. A., UV Spectroscopy of Planets, JPL TR No. 32-822 (December 15, 1965) (see Figure 20) 

(b) Ibid, see Figure 5 

(c) Ibid, see Figure 8 

<d) Ibid, see Table X 

(e) Ibid, see Table 

(f) Lawrence, G. M. , Journal of Quantitative Spectroscopy and Radiative Transfer, _5, No. 2, 359 (1965) 

(g) Dalby, F. W. , Handbuch der Physik, 27, 464 (1964) 

(h) Dufay, M. and Poulizae, M. C. , Memoires de La Societe Royales de Sciences de Liege, pp. 427 (c. 1966) 



species as OH, CN, N2 + , (>2 + , CO + , and CC>2 + layering could be invoked. 
Furthermore, for such molecules as CO, NO, and O 2 the atmospheric distri- 
bution question can not be answered at this time. Additionally, an im- 
portant factor in evaluation of twilight observation of molecular fluo- 
rescence involves the horizon screening height shown in Figure 42 and 
discussed in the next section. For example, if the initial absorption 
occurs at a wavelength of about 1400$, (i.e. , see the 4th positive 5-0 
transition for CO in Table 4) the resultant screening height is approx- 
imately 160 km. When this value is added to the h 2 -values of Figure 40, 
it is evident that the column count of solar-illuminated CO molecules is 
relatively limited. Accordingly, under certain conditions it may be 
advantageous to select absorption transitions which involve relatively 
long wavelengths and correspondingly lower screening heights. Addition- 
ally, it should be noted that unless the molecular constituent resides in 
an atmospheric layer configuration, the performance of twilight measure- 
ments may not be particularly advantageous. Alternatively, the perfor- 
mance of both day and twilight airglow experiments could provide data 
from which the existence of either the mixing or layering configuration 
could be ascertained. Finally, the role of quenching must be considered 
for relatively low altitude airglow phenomena. This is especially true 
in cases where forbidden transitions (i.e., long lifetimes) are involved. 

Many of the transitions specified in Table 4 can result in fluores- 
cence caused by the interaction of photoelectrons in the upper atmosphere 
of Mars. It should be noted that in the earth case for certain consti- 
tuents,' this phenomena represents the major emission factor so that its 
effect on Mars should also be recognized. However, this well known addi- 
tive source function is not stressed here ance it is a tractable problem 
so long as the appropriate data are available. 

On the basis of the foregoing discussions, and the data presented 
in Table 4, it is evident that the definition of specific measurements 
is an extremely complex matter. However, with appropriate effort, some 
selectivity for the ensemble of possibilities in Table 4 could be made in 
order to isolate the more promising cases. For the present purposes, this 
task is considered to lie beyond the scope of this presentation. 

In summary, then, it is felt that the data shown in Table 4 and the 
foregoing discussions indicate the potential of performing fluorescence 
experiments in the spectral region XX > 2000$, in order to identify a 
number of possibly important molecular constituents in the Martian atmos- 
phere and ionosphere. 

It is conjectured that various possibilities exist for the generation 
of chemiluminescent atomic and molecular emissions due to reactions between 
ambient chemical constituents in the Martian atmosphere. As an example 
of the former (atomic emissions), chemiluminescence may result .typically 
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from the debris species enumerated in Table 3 by mechanisms related to 
those which photochemically produce the nighttime sodium-D lines observed 
in the earth's atmosphere- In addition, twiiight and night airglow 
emissions at 6300/ 6364, 5199$, etc- may be produced from a variety of 
processes which include photodissociation, recombination, photoionization, 
electron and proton impact, etc. Concerning molecular system chemilum- 
inescence, it has been suggested that a significant portion of the terres-r 
trial night airglow is attributable to the production of produced 

by the reaction between NO and 0. This process as well as a similar one 
involving CO and 0 may occur in the Martian atmosphere. These latter two 
possibilities are discussed in a more detailed evaluation below where 
resultant order of magnitude signal intensity estimates are derived. 

Inn (Ref- 66) first suggested that the blue haze observed on Mars 
may be due to chemiluminescence from a reaction involving CO and 0 
which resulted in the production of excited C02*. The subsequent emission 
from this excited species resides in a band system in the spectral region 
XX 3200-6000$ with a peak intensity occurring at a wavelength of 4000A. 

The spectrum was photographed and classified by Mahan and Solo (Ref. 67) 
who reported that this system consisted of 250 bands with relatively 
little continuum background. Inn (Ref. 66) employed a photochemical 
equilibrium model (using local number densities for [CO] = [ 0 ] = 10^2 cm"3 
at an altitude of about 100 km) and the chemiluminescence efficiency 
data of Clyne and Thrush (Ref. 68) to estimate a total band system 
brightness of about 20 kR. On the basis of this value and the spectral 
distribution observed by Mahan and Solo (Ref. 67), it can be estimated 
that the signal brightness of the most intense band of the system could 
amount to about 1 kR. A signal intensity of this magnitude would be 
easily detectable under the appropriate twilight conditions (see Figures 
40 and 41) • 

With respect to the possibility of observing the N02 rf emissions in 
the twilight or night airglow, it is necessary to establish an upper 
limit value for NO in the Martian atmosphere. Marmo and Warneck (Ref. 69) 
have shown recently that under thermodynamic equilibrium an. N0/02 ratio 
greater than 1 obtains throughout the mixed region of the Martian atmosphere. 
Additionally, Belton and Hunten (Ref. 70) experimentally established an 
upper limit value for O 2 of about 20 cm-atm. or about 5.4 x 10^0 cm^ 
column in the Martian atmosphere. If this value is also applied to the 
case of NO, then a mixed atmosphere, at an altitude of about 100 km the 
local number density of NO, [NO] £ 10^ crti“^. Jonathan and Doherty (Ref. 

71) have demonstrated that the chemiluminescence efficiency for the N02* 
emissions is about 2000 times greater than that measured for the CO 2 * through- 
out the investigated spectral ranges. Thus, it is evident that under 
conditions similar to those specified by Inn (Ref. 66) (i.e., where [ 0 ] = 

1012 cm"^) the resultant chemiluminescent signal intensity may be as 
large as 4 x 10^ kR. However, it should be noted that this latter value 
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is predicated on the photon emission rate over the entire continuum which 
extends over a spectral interval of about 20008. Thus, on a unit spectral 
interval basis, the average chemiluminescence signal amounts to about 
2 kR/8, However, since the continuum peaks around 6200&, observations 
conducted at this spectral position can be expected to amount to about 
an order of magnitude greater or about 20 kR/A. Clearly, such signal 
magnitudes could be observed easily during the course of the proposed 
experiment. However, in reality the local NO number density at a Martian 
atmosphere altitude of 100 km is probably less than lO 11 - cra"^ so that a 
concommitant reduction in signal intensities would result. 

Finally, with respect to the two chemiluminescent reaction examples 
discussed above, it should be noted that the pertinent chemical consumption 
rates are such that the active constituent number densities are essentially 
conserved throughout the night. Under these minimum background intensity 
conditions, it is evident that a sensitive technique is available for 
positively observing and identifying these emissions and subsequently 
establishing meaningful upper limit values for. the constituents involved. 

In this section is described an experimental technique designed to 
detect, identify, and measure the altitude profiles of. minor constituent 
absorbing species in the lower atmosphere of Mars. For this purpose, 
twilight measurements are performed to measure the irradiances of both 
the zenith and horizon directions. ’ It will be demonstrated that these 
irradiances are sensitively perturbed by the presence of absorbing trace 
constituents in the lower atmosphere of Mars. This technique involves 
the following advantages: (a) the detection capability is enhanced by 
about a factor of 25 greater than that available with currently employed 
techniques; (b) the unambiguous identification of the constituent involved 
is possible owing to its absorbing characteristic for XX > 2000)8; and 
finally (c) the species altitude distribution can be ascertained on the 
basis of the geometric factors involved in the performance of cogent 
twilight experiments. 

As discussed previously, the Martian atmosphere is essentially 
transparent to solar zenith radiation for wavelengths above about 2000A 
as shown by the results of Figure 38 and confirmed by the recent observa- 
tions of Evans (Ref. 72). On this basis, it is possible to experimentally 
establish upper limit vertical column count values for any gas which ab- 
sorbs radiations of wavelengths greater than 2000.8, For example, employ- 
ing the measured absorption cross sections for ozone in this spectral 
region (Ref. 52), a corresponding upper limit vertical count value of 
about 1 x 1017 cm^ column can be established. Identification of the 
existence and establishment of the upper limit estimate of the content 
of this particular constituent has important implications on both aeronomic 
and biological processes occurring in the Martian atmosphere. 

Similarly, upper limit values can be established for other potentially 
important species which exhibit distinctive absorption features. The 



absorption curves for a number of such potentially important candidate 
species which qualify in this regard are presented in Figures 45, 46, and 
47 over the spectral region XX 2000-3000$. These data- Were extracted 
directly from Calvert and Pitts (Ref. 73). It should be noted that the 
curves pertinent to inorganic species are presented in Figure 45 which the 
absorption characteristics of several interesting organic materials, 
which may be related to biological interests, are shown in Figures 46 and 
47. 


With respect to the data of Figure 45, it should be noted that in 
biological processes, H 2 O 2 is required to oxidize both sulfides and 
certain carbon compounds to simple sulfates and carbonates. In addition 
Marmo and Wameck (Ref. 69) have demonstrated that SO 2 or H 2 S can parti- 
cipate as precursor sources for sulfur in the photochemical generation 
of sulfur containing amino acids in the Martian atmosphere. Marmo and 
Warneck (Ref. 69) have also shown that, in the Martian atmosphere, a thermo- 
dynamic relationship obtains for the NO/O 2 ratio, which can be further 
related to the content of NO 2 , N 2 O, and NOCi so that the presence of 
these constituents may involve both important aeronomic and biological 
implications. For inorganic gases, a number of additional candidate 
species have been enumerated previously in Table 4 where appropriate 
absorption (v n - 0) bands are shown to occur at X > 2000$ for several 
transitions. 

Absorbing gases (for XX > 2000$) which may be present in the lower 
Martian atmosphere may act as photosensitizers to allow solar photolysis 
of CO 2 .J HOH, organics, and other molecules which do not absorb the rela- 
tively intense solar flux in this spectral, region. For example, although 
CO 2 and water vapor do not absorb strongly in this region, the invocation 
of photosensitization by an absorbing species (i.e., S02> H 2 S) could 
result in the photolysis of these species and subsequent generation of . 
several simple .organic molecules which may have important implications 
on the exobiology of Mars. In fact, several compounds possibly produced 
by the above mechanism are also represented in the life-cycle. Character- 
istically, simple organic compounds often absorb strongly in the region 
XX 2000-3000$ as illustrated by the data of Figures 46 and 47 pertinent 
to a variety of formaldehydes, ketones, organic acids, hydrides, and 
esters. Calvert and Pitts (Ref. 73) present a number of photolytic 
mechanisms involving these and other organic materials, the products of 
which would undoubtedly be important to Martian biological investigations. 
The detection of these or other organic materials does not intrinsically 
assure a biological origin, although the presence of these constituents 
could represent and possibly be identified with key residue materials of 
biological processes. In this regard, it would be or particular value to 
review the experimental requirements and expected results of complemen- 
tary exobiological experiments aboard the Martian Lander. 

The salient features of an atmospheric absorption measurement experi- 
ment can be discussed both in terms of the geometry illustrated in Figure 
37 and the corresponding twilight* zenith and horizon radiances presented 
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Figure 47 Absorption spectra of selected organic molecules. 
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in^-Figures 40 and 41, respectively. The detection of minor constituents 
will be achieved by performing both zenith and horizon spectral scans 
(XX 2000-3000X) throughout twilight observation periods. It can be seen 
from Figure 37 that a significantly enhanced cm^ column count is involved 
along the path ACD. On this basis, the presence of minor constituents 
along this path would sensitively affect the observed zenith and horizon 
radiances as a function of both solar depression angle and wavelength . 

The application of appropriate analytic techniques should result in the 
extraction of the absorption characteristics of the lower atmosphere and, 
under certain conditions, specific -constituents may be identifiable. 

Finally, if the major absorption contribution over a given spectral inter- 
val can be ascribed to a specific constituent (i.e. NO, CO, OH, CN, etc.), 
that it might be possible to obtain the scale height or number density- 
altitude profile for that species. In any event, it is evident that the 
highly sensitive. absorption experiment described above can serve to estab- 
lish new identifications, upper limit content values and distributions of 
trace constituents resident in the lower Martian atmosphere. 

There appears to be little question regarding the occasional exis- 
tence of dust and/or haze layers in the lower Martian atmosphere (Ref. 74) 
Physical, photographic and polarometric measurements have indicated that 
several cloud forms exist including white clouds, blue clouds, "yellow 
veils", and "violet hazes". The white clouds probably involve a particle 
size range between 1 to 60 microns and are believed to consist of either 
ice or solid C02* It is believed that their formation is favored in low 
temperature regions since they have often been observed at the sunrise and 
sunset limbs of the planet. The blue clouds, characterized probably by 
a particle size range between 0. 1 to 1 micron have been compared by some 
to terrestrial noctilucent clouds, while others feel they represent an 
accumulation of micrometeorites suspended in the lower Martian atmosphere. 
The relatively rarely observed low-level yellow clouds probably consist 
of white desert dust blown about by Martian surface winds. In summary, 
it appears that not only is the presence of dust in the Martian atmosphere 
confirmed, but that the particle diameters appear to range from 0. 1 to 
60 microns. 

The general experimental configuration and the procedures applied 
herein should closely resemble those described in the previous section cm - 
- cerning detection of trace molecular constituents by absorption differen- 
tiation over the spectral region XX 2000-8000S. However, the conduct of 
the observations should not be confined to twilight since it certainly 
appears feasible to search for particulate scattering from dust layers 
under full solar illuminated daytime conditions as well as during times 
of Phobos illumination on the night side of the planet. 

For the ranges of wavelengths and particle diameters involved, it 
appears that both Rayleigh and Mie scattering contribute to the observed 
signals, so that performance of the spectral scanning procedure over the 
specified geometry should result in confirmation of the presence or' ab- 
scence of dust or haze layers in the lower Martian atmosphere. Additionally, 
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application of appropriate analytic procedures to representative spectral 
scans under the variety of conditions encountered during the experiment 
may produce data relevant to the size distribution, nature, and other 
physical characteristics of the particulate matter. Twilight observations 
could be employed for example, to locate the tops of layers. While this 
latter information is of great interest intrinsically, it could also be 
applied to the analysis of the other proposed twilight measurement exper- 
iments discussed previously. 

With respect to nighttime measurements, the white clouds (size 
distribution from 1 to 60 microns) appear to form in the vicinity of the 
planetary sunrise and sunset limbs, so that such clouds may exist through- 
out the entire dark side of the planet owing to the relatively low am- 
bient temperature conditions. This possibility could be verified simply 
by observing' the presence of Mie scattering (far more efficient than 
Rayleigh scattering) from the particles illuminated by Phobos during its 
orbit into 'and around the dark side of the planet. For this particular 
case, however, preliminary calculations indicate that the relatively low 
signal intensities involved would probably require that the spectral re- 
solution of the instrument be increased to about 100$. 

Finally, the employment of a polarization disc should be considered 
in order to acquire pertinent information to further characterize the 
size, distribution, and physical nature of the observed particles. How- 
ever, it should be noted that employment of such a polarization disc gen- 
erally involves signal intensity losses as well as imposing a short wave- 
length limitation on the observations. 

The discussions presented in the previous sections have identified 
the general operational requirements and constraints of the proposed in- 
strumentation system. Among the specific requirements are included (a) 
spectral scanning capability to cover both the ultraviolet and visible 
regions from \\ 2000-8000$; (b) a spectral resolution of approximately 
lA would be desirable; (c) a highly ruggedized instrument having minimum 
moving parts due to the nature of the mission; (d) sensitivity such that 
brightness values of between 1 to 10 R/$ can be observed over the ultra- 
violet spectral region; and finally (e) employment of a double-pass con- 
figuration in order to minimize significant instrumental scattering prob- 
lems which may be encountered in the adopted experimental mode. 

In this regard, under Contract NASS -947 2 the GCA Technology Division 
designed a double-pass UV scanning spectrometer for operation in a sat- 
ellite, environment for the purpose of measuring the content and distrib- 
ution of ozone in the earth atmosphere. The salient features of this 
latter instrument shown schematically in Figure 48 can be applied directly 
to the proposed experimental requirements enumerated above. The general 
features and pertinent operational aspects of the instrument are des- 
cribed elsewhere (Ref. 50). 
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Figure 48. Schematic diagram of proposed instrun^nt. 


Order of magnitude calculations indicate that the employment of a 
5 x 5 cm grating with 2400 lines /mm and a 50 cm focal length would result 
in the achievement of a 1A resolution with a slit width of about 25 microns 
(if the grating is set for an angle of incidence of 45 degrees as shown in 
the figure). These calculations also indicate that the resultant photon 
flux at the exit slit will be about 100 photons /second for an emission 
rate of ten R/A. Using these parameters , further analysis indicates that 
the entire spectral interval from 2000-8000& would be spread over a length 
of about 15 cm across the focal plane as shown in Figure 48. 

It should be noted that the optical components shown in Figure 48 
are schematically indfcated by lenses; however, in practice mirrors would 
be involved. There is no requirement for the optica-1 components to vary 
their positions during scanning operations so that the fabrication of a 
ruggedized version of this design should constitute a relatively straight- 
forward task. 

Much of the material generated under this study has been utilized in 
the preparation of a suggested experiment for the forthcoming Mars Lander 
program. 

"The Rocket Measurement of the 62 Micron 0 I Emission in the Upper 
Atmosphere." The elucidation of the various processes responsible for the 
thermal structure of the upper atmosphere of the so-called thermosphere 
(Ref. 75) is of importance for the understanding of the thermal structure 
of the planetary upper atmosphere. The thermal structure of the lower 
atmosphere is mainly controlled by the local absorption and emission of 
radiation, and the transport of thermal energy by conduction is relatively 
unimportant. However, Spitzer (Ref. 76) pointed out that, in the upper- 
regions of the atmosphere, the role of absorption and emission is very much 
decreased due to the low total particle density, and the conduction of 
thermal energy becomes important. This is due to the fact that the thermal 
conductivity of a gas is independent of the gas density. 

Bates (Refs. 77,78) has considered in detail the heat input by the 
absorption of solar ultraviolet radiation and its loss by the 62 micron 
emission by atomic oxygen. The thermal balance of the atmosphere is 
given by 

k x - -G + R, (29) 

dz 

where k is the thermal conductivity, 

R is the loss of heat by the radiation of 62 micron emission, and 
G is the local heating caused by solar uv radiation and is given 
by 

G(z) = 6 1 ("Oct n(Z) e * n ( Z > cr H (30 \ 
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where I(°°) is the solar flux at the top of the atmosphere, 
cr is the absorption cross section, 

n(Z) is the particle concentration of the absorbing constituents, 

H is the scale height, and 

€ is the heating efficiency^ factor or the fraction of the abs'orbed 
uv energy which is used in local heating. 


The information regarding solar VUV radiation flux and that of ab- 
sorption cross section. is fairly satisfactory. However, the heating 
efficiency factor is poorly known. McElroy, et al (Ref. 79) have dis- 
cussed the various processes responsible for the conversion of the ab- 
sorbed VUV radiation into local thermal energy. • Their discussion indicates 
that, from the consideration of the complexity of the processes, the 
theoretical determination of the heating efficiency factor is difficult. 
Hunt, et al (Ref. 80) have concluded that the estimated heating efficiency 
factor varies between 6 percent and 50 percent. It is proposed here 
that the simultaneous observation of the transmission of solar VUV rad- 
iation and that of 62p emission by rockets can yield reliable information 
regarding the heating efficiency factor. 


If we substitute Equation (30) in Equation (29) and integrate, we 

have 


k x 


dr 

dZ 


!(«>)€ [1 - exp (-n,(Z)cr H)) 


CP 

-/ 


R(Z) dZ.. , oi; 


At thermopciuse, tne tempera cure gradient is zero and; therefore, the 
thermal efficiency factor is given by 


€' 


oo 

/ 


R(Z)dZ 


i(«) ll - exp(-n(Z)cr H)] * 


(32) 


Therefore, it can be concluded that the measurement of 6 2(jl flux can give 
valuable information regarding the thermal structure of the atmosphere. 


The 0 I emiss.ion resulting 
structure levels of the ground ^ 
in Figure 49. 


from the transitions between the fine 
P state of the atomic oxygen are shown 
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kn(M)>> i^ A l-2 

» 5. 34 x 10 


-5 


or n(M) » 


5.34 x 10 


-5 


(36) 


TABLE 5 


State *eV 

A -1 
sec 

s 



exp ( 

— e /kT) 
n 



100°K 

500°K 

1000°K 

1500°K 

2000°K 


% 

0.028 

1.7xl0" 5 

1 

3.88xl0* 2 

5. 22x10 

7. 23xl0 -1 

8.02X10" 1 

8.52X10” 1 

8. 952 

3p i 

0.0.20 

8.9xl0~ 5 

3 

9.83xl0' 2 

6.31X10” 1 

7.95x10 

8.56X10" 1 

s^oxio -1 

9.23a 

% 

0.000 

— 7 ’ 

5 

... 

— 

— 

--- 


— 


Ig n exp(-6 n 

/kT) 

5.33 

7.41 

• 8.11 

'8.37 

8.52 

8.6 


The information 'regarding collisional energy transfer between the fine 
structure levels of the ground state atomic oxygen by process (33.) is not 
available. However , the information regarding the rate constant of simi- 
lar processes for Hg and alkali atoms (Mitchel and Zemansky) (Ref. 81) is 
known, and it varies between 5 x 10~10 cm3/sec for Na and 1.5 x 10 -i ^ for 
Hg. Assuming k ~ 10 cm^/sec, we find that the above condition is sat- 
isfied for 

n(M) » 5.3 x 10 7 . (37) 

Therefore, it can be concluded that the distribution of atomic oxygen be- 
tween the fine structure levels of its ground state will be in thermo- 
dynamic equilibrium with the ambient up to altitudes of about 250 km in 
earth's atmosphere; of course this may be correspondingly higher if the 
rate constant is larger than the value used above. 
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Figure 49 


The 3 P 0 and 3 P X levels of 0 I lie only about 0. 028eV and 0.020eV above 
the lowest 3 P 2 level, so that the excitation potentials involved are of 
the same order as the energy of the thermal motion. The transition pro- 
babilities of the two transitions 3 Pq -*■ 3 Pi is about 1.7 x 10 "5 sec" , 
and that of 3 Pi -»■ %2 is about 8.9 x 10" 3 sec"^. If the collisional 
excitation of 3 Pq and 3 P-. is made sufficiently rapid by the following 


0CP 2 ) + M • 


0 ( 3 P X or 3 P 0 ) + M, 


then the distribution of oxygen atoms in these levels will be 'the char- 
acteristic of the ambient temperature. Since the lifetime of 3 P^ level 


is shorter than that of 
iation. 


level, we shall mainly consider the 62{j, rad- 


The required condition for the Boltzmann distribution of 0 atoms be- 
tween 3 Pi and 3 ?2> which is characteristic of the ambient temperature, is 

k x n [0, 3 P 2 3 n.[M] exp (-e^kT) » n [0, ^ ] (34) 


. n [0, 3 P ] - — exp <-€ 7 /kT) n [o, 3 P ] , 
1 s 2 1 2 


o o 

where g^ and g 2 are the statistical weights of the two J P]_ and P 2 states, 
respectively. T is the ambient temperature and eq is the excitation energy 
of the 3 P^ level, which is equal to 0.02 eV. The relevant information is 
collected in Table 5. Using relations (34) and (35), the condition for 
the thermal equilibrium between 3 P^ and 3 P 3 can be written as 
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The loss of radiation in a planetary atmosphere by atomic oxygen 
depends on its optical thickness above the reference altitude, which in 
turn depends .on its absorption cross section. The .absorption cross sec- 
tion at the line center is given by 


_L 2 \ 112 £1 1 A 1 2 

a 2-I * 4X ' X > ' g 2 Av d 


wnere A is the wavelength of the line center and 


(38) 




2kT In 2 
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1 1/2 


1 

X I 


is the DoDoler ha If -width of the line; Therefore, 

,ll/2 


. “?-! * 4X 


me 


2XkT 


S 1 3 

— x X A 

8, A 


l->2 


= 9.4 x 10" 17 x t" 1 ^ 2 cm 2 ; 


. (39) 


and the optical thickness of the atmosphere is. given by 

’ CO * 

-/ a 2-l [1 " ex P (-e x /kT)] n <0, 3 P 2 ) dZ (401 

Z 

► 00 

't ~ 9 . 4 x 10~ 17 / » [1 - exp (-e /kT)] n (0, 3 P.) dZ 

^ sTt 1 2 

and for an isothermal atmosphere 

00 

U • £XP ( - € l /kT)] x g exp(-t 7k?) / n(0)dZ 
n n - 


<=9.4x 10‘ 17 x. 1 


00 

/ n( 


" = K x / n(0) dZ 


(41) 
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The calculated values of K' for different temperatures are given below: 


TABLE 6 . 


Temperature 

°K 


K 


100 


8.82 x 10 


■18 


500 


9.75 x 10 


1000 


•19 


3.76 x 10 


-19 


1500 


2.08 x 10 


•19 


2000 

1.36 x 10 


Temperature 

°K 3000 

K 7.63 x 10“ 20 


It can be concluded that, for the same optical thickness at the line 
center, larger numbers of atomic oxygen per cm^ column are required at 
higher temperatures. This is understandable because at higher temperatures 
the absorption cross section at the line center decreases and also the 
concentration of 0 ^P2 The rate of emission of radiation after 

neglecting the effects due to large optical thickness is given by 

R(0) s ^l-+2 n ^ x 1.6 x 10 erg/crrf* * sec 

- n x 1.6 x 10 ^ erg/cnf* * sec, (42) 

where T| is the fraction of total atomic oxygen concentration which are 


in the P state. 


g x expC-Gj/kT) 

^ ~ 8 0 + g x exp (-e^/kT) +g 0 exp(-£' 2 /kT) 
Hence, the net emission rate is given by 

3 

R(0) - R^, n(0) erg/cm • sec, 

where ^ x x r\ x 1.6 x 10* 12 

which depends on the temperature. Substituting from Table 5, 
culated Rj and its values are given below: 


(43) 

(44) 

(45) 

we have cal- 
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TABLE 7 


Rrp (Erg/atom • sec) 


Temperature °K 


100 

1 . 

.58 

X 

10- 19 

,500 

7. 

,27 

X 

lcf 19 

1000 

8. 

,38 

X 

ID ' 19 

1500 

.8. 

,75 

X 

IQ ' 19 

2000 

8. 

,'92 

X 

ID ’ 19 

3000 

9. 

,09 

X 

io “ 19 


The resultant flux of 62p. radiation in the earth's atmosphere can 
now -be calculated for a given standard model atmosphere. The relevant 
parameters for CIRA (1965) model are shown in Table 8. Using Equation 
(41) and the values of K given in Table 6 , the optical thickness of the 
atmosphere above different altitudes are calculated and shown in the last 
column of Table 8. It can be concluded from Table 8 that the optical 
.thickness for 62ji radiation is much less than 1 at altitudes above 150 
km. It should also be noted that the required condition (Equation 37) 
for the Maxwell-Bo It zmann distribution is satisfied at least 'up to 300 km 
and may be to higher altitudes. Therefore, the flux of 62p. emission rad- 
iation in a region between 150 km and 300 km can be calculated from 
Equations (44) and (45) and using the values of given in Table 7. 

The calculated rate of emission per unit volume is given in Table 9. The 
rate of emission per cra^ column seen by a detector pointing upward in a 
rocket is given in the last column of Table 9 and is obtained by multi- 
plying the volume emission rate by the atmospheric scale height. It 
should, however, be noted that the emission flux seen by a rocket borne 
detector pointed. in the horizontal direction is about 30 times that seen 
in the vertical direction. It can be concluded from Table 6 that-, be- 
tween 100 to 300 km altitude, the observed flux varies between about 10”& 
to 2.5 x 10“^ watt/cm^ for a detector pointed to the upward direction 
and between 3 x 10 to 7 x 10 watt/cm^ for a detector pointed to the 
horizontal direction. 

It is now appropriate to consider the feasibility of detecting the 
62p radiation by means of a rocket platform experiment. 

The detectability of the radiation by 0.1 at 62|i depends on celes- 
tail lackground at the top of the atmosphere due to stellar radiation and 
that originating in the solar system. .The infrared background at the top 
of the earth's atmosphere, due to the latter source, is due to the re- 
flection of solar radiation by planets and satellites and also due to the 
thermal emissions. A major component is also due to the infrared radiatior 
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TABLE 8 



300 


1755 


46.40 


TABLE 9 


Altitude 

Emission Rate 
Per Unit Volume 
Erg/cm 3 sec 

Emission Flux 
2 

Erg/cm sec 

100 

1.55 x 10“ 7 

9.98 x 10~ 2 

150 

1. 10 x 10~ 8 

2.53 x 10" 2 

200 

2.75 x 10" 9 

1.03 x 10~ 2 

250 

1.03 x 10~ 9 

4.97 x 10“ 3 

300 

4.52 x 10" 10 

2.55 x 10~ 3 
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from zodiacal light and gegenschein. The information regarding the infra- 
red background at the top of the earth's atmosphere is not available in 
the 62|r region. However., the total infrared irradiance in the 0.8 to 3p 
region is as follows (Handbook of Geophysics and Space Environ" 10 '"*- 0 1Qfi ^ 


Source 

Stellar 


Irradiance 

W/cm 2 



Sun 

1. 19 

X 

10‘ 

-l 

Moon (Full) 

2.6 

X 

10' 

-7 

Mercury (Elongation) 

1.2 

X 

10* 

*11 

Venus (Elongation) 

1.8 

X 

io* 

-10 

Mars (Quadrature) 

3.7 

X 

10' 

•12 

Jupiter (Quadrature) 

8.2 

X 

10~ 

•12 

Saturn (Quadrature) 

8.4 

X 

io" 

•13 

Uranus (Apposition) 

1 

X 

10" 

■15 

Neptune (Apposition) 

S . 

X 

io7 

■15 

Zodiacal and Gegenschein 


in 


(between Elongation of 50 

1 6 ' 

• iU 



It can be concluded from the above table that the infrared irrad- 
iance between 0.8 to 3ji at the top of the earth's atmosphere and in the 
absence of solar of lunar illumination is about • 10”^ W/cm^. If this 
radiation is assumed to be characteristic of effective black body temp- 
erature of a few thousand degrees (the effective black body temperature 
of .the sun and other stellar sources) then at longer wavelengths, the 
irradiance should decrease as X~^*. Therefore, the extrapolated irrad- 
iance at the top of the earth's atmosphere is about 10" ^5 w/cm^ around 
62fi 0 I emission is several orders of magnitude above the natural back- 
ground due to extra terrestrial sources. 


We shall now consider the detectability of the 62p. 0 I emission by 
an appropriate rocket-borne detector for this spectral region. A suit- 
able detector in this region is the. low- temp erature gallium-doped ger- 
manium bolometer system which operates a liquid helium temperature (Ref. 
82). The general characteristics of this detector have been discussed 
elsewhere (Ref. 83) and are not repeated here. Specifically, the total 
flux Fp received by the detector is approximately given by 
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where: 


D = ?! / n(0) dZ 

z 




D.= element area 

2 

F = emission flux (in W/cm ) 

SI = effective solid angle of the instrument 
K = instrument constant 


The N.E.P. of the above types of detectors is about lO"^ watts 
(Ref. 82), and it can be concluded ' that -the atmospheric emission at 62^i 
of 0 X can be successfully measured by available detectors since the 
emitted flux is between 6 x 10“ W/cm^ to 10“ 10 W/cm . 

"The Deposition of' Meteoric Debris Constituents into the Upper Atmos- 
pheres of Earth, Venus and Mars." Notwithstanding the lack of quantita- 
tive data on the inf lux. rate and mass distribution of meteoric debris 
sufficient data are available to make it evident that the major contribu- 
tion to the accretion of this debris by earth is due to particles with 
masses less than about 10“ 7 grams (which corresponds to radii less than 
50{j.) . (Refs. 84-86). Prior studies (Refs. 87-95) to determine flux from 
rocket collections have assumed that no mass loss by ablation occurs 
when the radius of the particles is less than 50p. A series of studies 
on the interaction of micrometeoroids with, the earth atmosphere (Refs. 96-98) 
have made it -evident that the above assumption is not true and that par- 
ticles with radii less than 50[i do indeed ablate so as to result in a 
deposition of micrometeoroid constituents into the upper atmosphere of 
earth. Specifically; Whipple (Ref. 97) demonstrated that spherical par- 
ticles of density 3 grams/cm^ entering vertically with velocities of 12; 

15; 30; and 70 km/sec ablate for radii less than 39; 19, 3; and 0.3p, 
respectively. Clearly; the deposition of this debris into the upper 
atmosphere of planets. can have important aeronomic effects as well as 
contributing to ionospheric phenomena; noctilucent clouds; haze layers; 
airglow phenomena; etc. It is the purpose of this study to obtain estim- 
ates of the mass' deposition rates of ablated atomic matter within the 
atmospheres of earth; Venus; and Mars. In accordance with the above 
discussion; this study is limited to meteoric debris particles with masses 
less than 10“7 grams; i.e. radii less than 50p. 

Allen; et al. (Ref. 99) have shown that shielding and fragmentation 
do not obtain for all r-values less than 50|i; whereas Levin (Ref. 100) 
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has shown that heat storage effects can be neglected for such particles 
so long as the entering velocity is less than approximately 60 km/sec. 
Accordingly, for this preliminary study it was felt that the classical 
theory of meteor ablation (Refs-. 96-98, 101-105) could apply in order to 
obtain the desired estimates. The appropriate techniques are first dev- 
eloped and applied to the atmosphere of earth. These findings are then 
extrapolated and applied to appropriate model atmospheres for Venus and 
Mars in order to obtain corresponding ablation-mass profiles. 

According to the classical theory for meteoric ablation (Ref. 106) 
once the ablation temperature has been achieved the instantaneous mass, m, 
of the ablated meteoroid at altitude, z. 


yFv 

8‘pr^cosX. 




where: = the initial mass 

7 = the drag coefficient (a value of 1. 1 is assumed here 

after Reference 107) 

2 

F = the shape factor, [defined as cr/jtr where cr is the cross 
sectional area of ’ the meteoroid arid r is the radius of 
a spherical meteoroid (F = 1 for a sphere)] 

v = the velocity 

11 2 2 

£ = the energy of ablation (a value of 1. 1 x 10 cm /sec 

is assumed here after Reference 107) 

p = the meteoroid density 

r = the initial radius (of a spherical meteoroid) 

X - the zenith angle of the radiant 
=■• the atmospheric density 


For ready solution of Equation (46) the atmosphere is considered to 
be composed of a series of concentric layers wherein each zone is suffic- 
iently small compared to the scale height, H. so that for computation 
purposes an average scale height H, obtains. Under these conditions. 
Equation (46) can be approximated by 


yFv 2 H 


8£pr cosX 


( p a - p 0 ) 


(47) 


which obtains for each individual layer. 
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The corresponding ablation rate, n, or atoms ablated per unit time 
for an individual particle is then given by: 


n 


1. dm 

II dt 


3 n o p a 7Fv3 

8 £P r o 


1 - 


yFv 2 ^ 


8£pr cos X ^ p a”°o^ 


— 1 3 

04 

atoms/sec 


or 


3n’p XPv 2 

o r a 

8gpr Q cos X 


1 


&ip r 0 cos X ^a~ p o^ 


(49) 


atoms deposited/unit vertical path 


The the above equations the subscripts Q refer to the upper boundary 
of a given zone and n 0 ' is the number of unablated atoms remaining in the 
particle at the upper boundary of a given zone and p. is the molecular mass 
(grams /molecule) of the meteoric material. 

The present study is concerned with the chemical interaction of the 
atomic debris with atmospheric constituents so that the role of the sur- 
viving particles is not considered further. However > it is clear that 
these studies can be applied to the latter problem for ‘the purposes of 
relating incoming flux values with empirical data obtained by satellites, 
rockets, balloons and surface measurements. In order to apply these equa- 
tions to the present problem, a number of meteoric physical parameters 
must be established. For this purpose the following comments apply. 

For the present purpose the particle influx-particle mass relationship 
(Ref. 84) shown in Figure 50 is selected and considered to represent' an upper 
limit value for the particle influx rate. The indicated distribution has 
not been firmly established so that the present results should be viewed 
with appropriate caution. In any case, this distribution makes it evi- 
ent that the major contribution to mass deposition in the earth atmos- 
phere is indeed^ due to particles with masses less than approximately 10 7 
and 10 grams, i.e. less than 50p.. Since there is considerable disagree- 
ment in the literature on the correct particle influx-particle mass rela- 
tionship, it is appropriate that the techniques illustrated in the present 
study be applied to other distributions so that some upper and lower limit 
values can be acquired. 
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Concerning the assumed shape of the micrometeoroid it is appropriate 
to assume spherical bodies since these result in a considerable simplifi- 
cation of the calculations involved. However,, it is recognized that photo- 
graphs of these particles collected during rocket investigations (Ref. 87) 
show that the micrometeoroids can assume a variety of shapes which appear 
to be rather irregular. 

For the present purposes the particle density-particle radius rela- 
tionship shown in Figure 51 was adopted on the basis of fitting and extrapo 
ating a smooth curve to the several derived densities and estimates of den- 
sities appearing in the literature (Refs. 87, 107, 108). In this regard 
it will be later demonstrated that most of the mass deposition is due to 
the ablation of micrometeoroids with r-values between 0..20 and 4. Op, so that 
particle densities of about 2.0 to 3.5 grams/cm^ obtain. However, -here 
again it is appropriate to caution that this relationship has not been 
firmly established. 

A representative value of 23 grams /molecule is taken for the atomic 
weight of micrometeoroid particles since this is the average value for 
stony meteorites (Ref. 102) which constitutes 90 percent of all incident 
micrometeoroids . 

Perhaps the most important single parameter to establish is that of 
the velocity distribution of the incident micrometeoroid particles. The 
distribution adopted for this study leaned heavily on the data relevant to 
twelfth visual magnitude radar meteors established by Eshleman, et al (Ref. 
109). This procedure was adopted since there are no pertinent data on the 
size distributions considered here and the data cited represent the faint- 
est meteors for which velocity histograms have been established. The vel- 
ocity histograms reported by Eshleman, et al (Ref. 109) included micro- 
meteoroids c£ interstellar origin, i.e. velocities 'greater than -72 km/sec. 
However, since this component amounts to only between 1 and 2 percent of 
the total incoming flux (Ref. 110) it has been neglected in the present 
study. - Furthermore, the published velocities were slightly redistributed 
so as to include velocities between 11-35 km/ sec, which were absent from 
the published data. In -any case, the velocity histogram shown In Figure- 
52 obtains for the present study and is considered to be a reasonable 
representation of the distribution involved. Corresponding velocity his- 
tograms for the planets Venus and Mars have been derived on the basis of 
data shown in Table 10. 


TABLE 10 

Meteoroid Velocity Ranges for Earth and Proximate Planets 


Planet 

Mean Orbital 
Velocity 

v (sun) 

esc 

v esc (planet) 

v . . 
minimum 

V maximum 

Venus 

35 

49 

10 

10 

85 

Earth 

30 

42 

11 

11 

73 

Mars 

24 

34 

5 

5 

58 
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PERCENTAGE OF ALL V - VALUES 



10 20 30 40 50 60 70 


PARTICLE IMPACT VELOCITY- V^£ km/ s**): 

igure 52 Adopted distribution for percentage or all incident v-values 
vs. particle impact velocity for planet earth. \X7 




The resultant histograms are graphically shown in Figures 53 and 54 
for Venus and Mars respectively. 

It has been established (Refs- 96-98,101-105) that a critical mass/ 
surface area ratio (m/s) obtains wherein incoming micrometeoroids with 
m/s -values less than critical will survive passage through the atmosphere 
(with no ablation) owing to the fact that they can not attain melting temp- 
erature, i.e. about 1800°K for stony meteors (Ref. 98). Specifically, 
Whipple (Ref. 97) has shown that for a non-iso thermal atmosphere 


3eyT. 
m 1 

■? 

> 



rH 

<=a(T 4 -T 4 ) 

s (v+q) cosX 

W i 

L Av 3 

^ CO — * 


q/ (v+q) 


(50) 


Where: e - natural logarithm, 2.718 

= temperature at base of region for which temperature 
. gradient is v = dT/dz 

= atmospheric density at base of region 

q = pg/R, p « mean molecular weight of air, gms/mole, 
g = acceleration due to gravity, R = universal gas 
constant, 8.31 x 10^ ergs/mole-deg 

e - emissivity 

a - Stefan-Boltzman constant, 5.67 x 10 erg 

T = instantaneous temperature of micrometeoroid 

T q = exoat mo spheric temperature of micrometeoroid 

For a spherical particle ~ ^ (51) 

From Equations (50) and (51), rp-values were calculated for the re- 
spective velocity ranges at the three planets using the parametric values 
in Table 11. 

Using the relationship between radius and density shown in Figure 52, 
the data in Figure 55 were obtained to indicate the minimum threshold 
radii for melting as a function of initial particle velocity. Specifically, 
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PERCENTAGE OF ALL V„ - VALUES 


0IC226 * SOS 



Figure 53 Adopted distribution for percentage of all incident v-values 
vs . particle impact velocity for planet Venus. 
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22 


0JC226 *40! 



PARTICLE IMPACT VELOCITY -V^km/sec) 

Figure 54. .Adopted distribution for percentage of all incident v-values 
■ vs. particle impact velocity for planet Mars. 
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TABLE U 



Earth 

Venus 

Mars 

Bef erence 

7 

1.1 

1.1 

1. 1 

(57) 

cosX 

0.8 

0.8 

0.8 

(61) 

T 

o 

300°K 

' 300°K 

300°K 

(47) 

T 

1800°K 

1800°K 

1800°K 

(52). 

A 

1 

1 

1 

(54) 

£ 

30 

44 

38.3 

Model Atmosphere 

q 

3.43xl0“ 4 
deg /cm 

4. 48xl0 -4 
deg/cm 

1. 65xl0" 4 
deg/cm 

Model Atmosphere 

ref. 

alt. 

93km 

110km 

110km 

Model Atmosphere 

h 

181°K 

210°K 

172°K 

Model Atmosphere 

p i 

1.76xlO~ 9 

gm/cm^ 

1.22x10" 10 
gm/cm^ 

0. 932x10" 10 
gm/cm^ 

Model Atmosphere 

V 

4. 58xl0 -5 
deg/cm 

. 9xl0" 5 

deg/ cm 

3.73xl0" 5 

deg/cm 

Model Atmosphere 


•k 

The Model Atmosphere for Venus follows McElroy (Ref. 112 ) , and for Mars., 
House, et al (Ref. 113). 


121 



PARTICLE RADIUS (microns) 





then, the area to the left hand side of a curve defines the radius -vel- 
ocity domain of micrometeors, i.e., particles that do ablate upon entry. 

It can be shown (Ref. 97) that for an initial m/s value and a given 
velocity, one can compute the corresponding altitude at which significant 
ablation commences, i.e., the temperature of melting is first achieved. 

The appropriate calculations have been made for an exoatmospheric micro- 
meteoroid temperature of 300°K (Ref. 102) and a temperature of melting- 
ablation of 1800°K previously established for stony micrometeoroids (Ref. 
98). The resultant data for the planets Venus, earth and Mars are shown 
in Figure 56. 

For the calculations to be performed on the interaction of micro- 
meteoroids with the atmospheres of earth, Venus and Mars the micrometeoroids 
mass distribution shown in Figure 50 is assumed to apply for each case. 

This is justified on the basis that no relevant data exist for the other 
planets. On the other hand, it is evident that the particle influx rate 
can be expected to vary for the three planets; however, in the present 
approach the debris deposition rate altitude profiles are normalized to 
the total incident ■ flux so that the resultant values can be readily adjusted 
in a linear manner on the basis of specific values adopted* Once this is 
established, it is a straight-forward task to demonstrate that more than 
99 percent of the incident mass ablated into a planetary atmosphere re- 
sults from micrometeoroids whose radii are less than only 4. Op.. To dem- 
onstrate this the data in Figures 50 and 52 have been convolved to deter- 
mine the percent of incident mass-incident particle radius relationships 
as shown in the histogram of Figure 57. At this point it is appropriate 
to summarize the following factors: (1) more than 99 percent of the inci- 
dent mass is due to raicrometeoroids with r-values less than 4. Op, (2) for 
a given r-value only those raicrometeoroids entering at velocities above 
the threshold value will significantly ablate (see Figure 55), '(3) for 
all incoming particles with m/s-values greder than critical -there is an 
altitude at which ablation commences which is a function of the particle 
impact velocity (see Figure 56). 

Figures 52, 53, and 54 represent plots of the percentage of incoming 
Vw-values against particle impact velocity for the planets earth, Venus 
and Mars respectively. It is a straight-forward task to compute the 
fraction of mass that will ablate for any given velocity increment in 
order to derive some idea of the fraction of material ablating in a given 
atmosphere. The relevant histograms are shown in Figures 58, 59 and 60 
for the planets earth, Venus, and Mars. In these histograms the .shaded 
portions represent the ablated fraction in a given velocity increment. 

Thus, the integrated shaded area in each case gives a measure of the per- 
cent of incident meteoroid mass that ablates within a given planetary 
atmosphere. These percentages are 86, 87 and 58 for earth, “Venus and 
Mars, respectively. An analogous calculation has been performed for the 
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ALTITUDE (km) 


OIC226 -! 



PARTICLE IMPACT VELOCITY - V ( km/sec ) 

_ — 

Figure 56‘. Threshold altitudes vs. particle impact velocities where' 

incoming particles commence melting at a temperature of 180£ 
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INCIDENT PARTICLE RADIUS Imitiwii) 

Figure 57 < Plot of percent of incident me** contained within, incident 

particle radii increment# between 0*2 and 0*8 micron# ■ fhie 
plot 1# derived on the Wi* of the di#cr ifcjt ion #hown in 
Figure 50* 
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f igure 59 . Plot ot velocity distribution »s shown Ip, Figure 53 wherein 
the percent of abUtcd vtniil in a giver velocity incre- 
ment la *hwn by the shaded regions. 



PARTICLE IMPACT VELOCITY -V ( lim /sec ) 


PERCENTAGE OF ALL V - VALUES 




planet earth to derive a relationship between the percent of deposited 
mass and the incident particle radius* In this case for each Individual 
particle radius it was assumed that the velocity distribution shown in 
Figure 52 obtains. The resultant histogram is shown in Figure 61 by the 
shaded portion which is superimposed on the data shown in Figure 57 for 
convenient reference* This display makes evident the fact that for micro - 
meteoroid particles with revalues less than 0.2u only an insignificant 
fraction of incident material is ablated into the upper atmosphere* This 
behavior of course , serves to limit the contributing radii to values be- 
tween 0* 2-4*0^* 


It is now appropriate to calculate the ablation rate-altitude profiles 
for the planets earth, Venus, and Mars* For this purpose it is convenient 
to rearrange Equation <49 ) as follows: 


3F7v^Tf 

cosX 



As 


Fyv 2 ? 


B^pr^ cosX 



(52) 


Where: v - the average velocity for a velocity increment 

r Q - the average radius for a radius increment Ar 

f = (Am/Z Am) ( Av/Z Av) is the convolved fraction of (a) the 
ratio of the incident mass contained In the radius incre- 
ment, Ar, to the total incident mass and (b) the ratio 
of the histogram percentage of velocities in the appro- 
priate velocity increment -- above the threshold velocity 
for the corresponding value of r — to all meteoroid 
velocities (100 percent)* 


For computational purpose to each velocity increment an average v -value 
was assigned and a series of mass deposition-altitude profiles were then 
computed for a given radius increment, average value radius increment, Y q * 

The final curves of altitude vs. normalized ablation rate are shown 
in Figure 62 for the three planets earth.* Venus and Mars* It is evident 
that these new results are predicated on the fact that the major deposition 
process is due to micrometeoroid particles with r-values less than about 
4*<Vu Hence, previous investigators (Refs. 101*114) employed earlier 
models of meteor influx rates which more heavily weighted contributions 
from larger particles so that the calculated peak altitudes were consider- 
ably lower than those shown by the solid curve (earth values) In Figure 62* 
On the basis of these published data (Refs- 101,114) Marino and Brown (Ref. 
115) arbitrarily adopted a peak deposition altitude of 105 km for their 
investigation on the role of interplanetary debris in planetary atmospheres* 
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ALTITUDE (km) 




Finally, these new data can find direct applicability to the mass 
spectrometric observations of metallic ions (presumably due to meteoric 
sources) recently reported by Narcissi, et al (Ref. 116), it should be 
recognized that the calculations performed in the present study are rather 
speculative and should be used with appropriate caution. However, they do 
serve to illustrate the general nature of the deposition of meteoric debris 
in planetary atmospheres so that they may be applied with some degree of 
confidence to some of the aeronomic problems involved in the interpreta- 
tion of observational data in the upper atmosphere of earth, Venus, and 
Mars. 


"Solar EUV Induced Fluorescent Radiation in the Atmospheres of Mars 
and Venus." Using photoionization production rates for excited states of 
C02"^ estimates have been obtained for fluorescent radiation from the 
atmospheres of Mars and Venus. The- most prominent spectral feature ex- 
pected is the band emission at 2882 and 2896.8, which may have an intensity 
as high as 8' kilorayleighs in the atmosphere of Mars and 35 kilorayleighs 
in the atmosphere of Venus. 

Using model atmospheres based on the pure CO 2 atmospheres of McElroy 
(Ref. 117) photoionization production rates have been determined for the 
X ^Hg, A b and C + states of C0 2 + . Ipn production by photo- 

electrons has been estimated by assuming that electrons with energy be- 
tween 34 and 68 eV produce one ion pair and electrons with higher energy 
produce one ion pair per 34 eV. The total ionization was divided among 
the four states of C0 2 +- using the experimental data of Bahr, et al. (Ref. 
118) for wavelengths greater than 6008. The division at shorter wave- 
lengths and that for electron induced ionization were determined by assum- 
ing the relative values at 6008 remain constant at shorter wavelengths. 

The total optical and photoionization cross sections used are those given 
by Henry and McElroy (Ref. 119). 

The accompanying drawings show the total ionization production rates 
and the integrated production rate above a given altitude for solar zenith 
angles of 0, 60, and 70 degress (See Figures 63-70 inclusive). 

In the absence of any deactivation process, each ion produced in the 
A 2]]^ state will decay to the X state by band emission at wave-lengths 
greater than 30008. A large fraction of the ions produced in the B 
state will decay to the ground state by band emission at 2882 and 2986S. 
Emission in these two bands will probably be the most prominent dayglow 
feature in the atmospheres of Mars and Venus. The wavelengths and inten- 
sities of emission from the C state are uncertain. 

Table 12 summarizes ion production computations. Table 13 presents 
total possible emission rates for each of the excited states, assuming 
no -deactivation or cascading. 
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TABLE 12 


Total Production Rate of CO^ States., Ions-cm 


Solar Zenith 






Venus Angle, Degrees 

X 

A 

B 

C 

Total 

0 

4. 55x10 10 

2.03 

3.53 

1.05 

11.16 

30 

3.94 

1.76 

3.06 

0.908 

9.67 

50 

2.92 

1.31 

2.27 

0.67 

7.18 

60 

2.28 

1.02 

■ 1.76 

0.52 

5.58 

70 

1.56 

0. 70 

1.21 

0.36 

3.82 

80 

0.79 

0.35 . 

0.61 

0.18 

1.94 


Mars 0 

10. 25xl0 9 

4.58 

7.95 

2.36 

25.15 

30 

8.88 

3.97 

6.88 

2.05 

21.78 

50 

6.59 

2.95 

5.11 

1.52 

16. 16 

60 

5.13 

2.29 

3.97 

1.18 

12.57 

70 

3.51 

1.57 

2.72 

0.81 

8.60 

80 

1.78 

0.80 

1.38 

0 

0.41 

4-37 


NOTE: Approximately 14 percent of the X il^ state ions are produced 

and approximately 25 percent of ions in the other states are 


produced in this x>?ay. 
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TABLE 13 



Emission Rate 

from CO* 

States, kiloRayleighs 





STATE 


Venus 

Solar Zenith 
Angle, Degrees 

A 

B 

C 


0 

20.3 

35.3 

0.5 


30 

17.6 

30.6 

9.-1 


50 

13.1 

22.7 

6.7 


60 

10.2 

17.6 

5.2 


70 

7.0 

12.1 

3.6 


80 

3.5 

6.1 

1.8 


Mars 0 

4.58 

7.95 

2.36 

30 

3.97 

6.88 

2.05 

50 

2.95 

5.11 

1.52 

60 

2.29 

3.97 

1.18 

70 

1.57 

2.72 . 

0.81 

80 

0.80 

1.38 

0.41 


NOTE: Neglects any cascading, de -activation 
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"Spectrography of the Lunar Atmosphere.” The proposed experimental 
investigation involves the employment of a compact, light weight, auto- 
matically operating spectrograph on the lunar surface to photographically 
record signature spectra of resonance and/or fluorescence scattering 
from- solar illuminated lunar atmospheric species (or LEM contaminant species). 
The relatively high sensitivity achievable by the proposed instrumentation 
(ppm or less) and the extensive spectral region included in the .measurement 
program (XX 500 to 4500$) result in the conduct of a survey experiment on 
establishing the presence or absence of specific lunar atmospheric species . 

The proposed experiment is ideally suited for operation on the lunar . 
surface. First, the spectrograph film format results in simultaneous inte- 
gration over each individual spectral resolution so that the required 
broad spectral range can be covered at full sensitivity. In addition, 
since the proposed instrumentation encompasses a self-contained operational 
entity, power,' telemetry,- data sharing, interference, etc. characteristics 
need be considered in the conduct of the experiment. Upon comparison of 
the overall capabilities of the experiment performed at a remote obser- 
vation platform to the suggested lunar surface program, the following ad- 
vantages are apparent: (a) the experimental conditions remain essentially 
constant throughout the viewing time, i.e., the incident solar flux, the 
ambient background, column of specific species (signal source), the view- 
ing position and angle, etc. ; (b) species number densities are maximum 
at the surface so that the resultant cm^-column count is optimized. This 
feature is especially evident in the detection of relatively heavy con- 
stituents, such as krypton, xenon, argon, etc. (which have relatively 
low scale heights from a remote platform); (c) the required astronaut 
activity is simple, minimum, and vital; (d) the delection of the surface- 
based viewing .geometry is .extremely flexible; (e) the instrumentation can 
be located in the vicinity of the LEM vehicle so that an opportunity 
exists for detecting fuel contaminants and measuring their inventories 
as a function of time. 

The data acquired from the indicated experiments can be employed in 
the identification of dominant lunar atmospheric source and sink functions 
as well as the pertinent mechanisms involved. Additionally, the spectro- 
graphic data can be correlated with other Apollo mission measurement pro- 
grams results, i.e., solar wind indicators, seismology, meteoric impact, etc. 

Several observational experimental programs have been conducted pre- 
viously to ascertain the existence, constituency, and extent of a lunar 
atmosphere. A concise summary of these results involving a broad variety 
of experimental techniques has been presented by Sytinskaya (Ref. 120). 

Highly sensitive experiments were performed by Lyot and Do Ilf us (Ref. 121) 
and Dollfus (Ref. 122) who were searching for the effects of Rayleigh 
scattering from a lunar atmosphere located just beyond the horns of a 
crescent moon. They employed a corona graph to minimize instrumental 
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scattering of disk light and took advantage of the polarization character- 
istics of Rayleigh scattered light. Their results indicated that the 
total surface number density of the lunar atmosphere must be less than 
10“9 that of the sea level terrestrial atmosphere (i.e. less than 3 x 10 
era" 3 ). 

On the basis of the results of a radio refraction experiment involv- 
ing lunar occultation of the Crab Nebula, Elsmore (Ref. 123) deduced that 
a surface electron density of about 10 cm"^ in excess of that of the 
surrounding interplanetary medium obtains for the lunar atmosphere. On 
the basis of an extrapolation of these data ; he suggested a new upper 
limit value for neutral species at the lunar surface of only about 2 x 
10**^ earth atmospheres. 

Her zb erg (Ref. 124) suggested that the observation of characteristic 
fluorescences for specific species offered a sensitive species detection 
possibility. Subsequently, Kuiper (Ref. 125) established that lunar 
atmospheric SO2 amounted to less than 7 x 10"^-® earth atmospheres on the 
basis of the absence of absorption bands on observed lunar disk light. 
Recently, Marmo and Engelman (Ref. 126) defined a manned orbital experi- 
ment involving the observation of resonance radiations from solar illum- 
inated neutral or ionic constituents in the lunar atmosphere. It' was 
demonstrated' both that a number of potentially important resonance lines 
occurred in the visible spectral region and that a meaningful experiment 
could be performed at an earth -based observatory. 

Finally; Marmo and Manring (Ref. 127) utilized the Sacramento Peak 
coronagraph to image, occult, and track the lunar disk in an attempt to 
observe atmospheric contents of Nal, Cal, and Gall along the lunar lim|j>. 
The results indicated upper limit values for these species of about 10 / , 
cm^-column. 

For the purposes of the proposed experiment, the state of present 
development over the entire field of activity is sufficient so that no 
new technology is envisioned. For example, the required astronaut activ- 
ity is s tr a ightf orward though vital in nature, thus requiring minimum 
effort and instruction. Additionally, a prototype of the suggested 
instrumentation has been employed in a number and variety of high altitude 
rocket experimental programs designed to measure vacuum ultraviolet (VUV) 
air glow in the earth atmosphere. For the present application, the only 
requirements involve simple instrumental size reduction, operational com- 
patibility with the lunar environment, and the design of a removable 
film casette for subsequent return to earth. 

The Apollo mission experiment discussed herein involves the employ- 
ment of a simple, light weight, automatically operated, dual spectro- 
graph designed to perform observations in two spectral regions, i. e. • 

XX 2500-4500& and XX 500-2500&. The experiment is to be performed on the 
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lunar surface in the protective shade of the LEM vehicle. The astronaut 
will initiate the pre-programmed experimental sequence of tasks which th« 
instrument is designed to perform automatically. Subsequent to auto- 
matic shutdown, the astronaut will remove the capsule containing the rol; 
film spectroscopic data for return to earth for development, data reduc- 
tion, and analysis and interpretation. This initial experiment should be 
considered in terms of a general survey in that an extensive spectral 
region will be covered wherein the presence or absence of signals owing 
to resonance scattering from solar illuminated species in the lunar at- 
mosphere can be detected with relatively high sensitivity. The recommended 
experimental technique will yield data which can represent meaningful 
inputs to the definition of relevant mechanisms producing the measured 
species number densities. The proposed experiment will be performed with 
existing state-of-the-art instrumentation and, minimum, although vital, 
astronaut activity. The experimental objectives and results can be applied 
and correlated with other experiments performed under the Apollo mission - 
program. Finally, on the basis of the derived experimental results, the 
location of a permanent astronomical airglow observing site on the lunar 
surface may be indicated which would operate essentially continuously and 
in direct communication with the earth. 

A meaningful experiment designed to determine the constituency and 
extent of the lunar atmosphere must operate within the constraints of 
the accepted upper limit values of 10® cm“3 and 1C)3 cm - 3 for total sur- 
face number- densities for neutral and ionic species, respectively (Refs. 
120-123). Within these limits, several theoretical investigations have 
been performed, the results of which are described in the literature 
(Refs. 126, 128-131) so that a detailed discussion is not included here- 
in. Suffice it to note that a number of representative possible lunar 
atmospheric constituents have been selected for the purposes of illus- 
trating the present experimental capability and its applicability to the 
present state of knowledge. This species tabulation is presented in 
column 1 of Table 14 with corresponding relevant resonance wavelengths. 
Present theory does not allow one to estimate the inventories of the 
individual species which may be resident in the lunar atmosphere which 
represents the general rationale for the importance of the proposed gen- 
eral survey experiment. In lieu of this specific knowledge, the expected 
signal brightness values in Rayleighs (where 1 Rayleigh = 10^ photons /cm^ 
sec emitted in 4ft steradians), are tabulated in column 3 of Table 14 for 
surface number densities, n Q , of only 1 cm" 3 for the individual consti- 
tuents. This sensitivity is expressed in an alternative fashion in 
column 4 of Table 14 in terms of minimum detectable n -values for an 
assumed detectability of 1.0 Rayleigh. Finally, the last column of 
Table 14 lists the major possible source functions of individual consti- 
tuents considered as extracted from the literature (Refs 122-126). The 
above tabulation has been separated into two spectral categories, i.e. 

2500-4500&, and XX 500-2500&, since the proposed investigation will 
be performed using two spectrograph systems consistent with this spectral 


145 



TABLE 14 


SELECTED 

LUNAR ATMOSPHERIC 

SPECIES, POSSIBILITIES AND ASSOCIATED SPECTROGRAPHIC DETECTION CAPABILITIES 



Expected Signal Brightness 

Minimum 

Major Sources of Con- 


X(8) 

($, Rayleighs) for Surface 

Detectable n 0 - 

stituent in Lunar 
Atmosphere 

Species 

Number Density (n 0 ,cm"^)-l 

Value for o=lR 



XX 2500-4500$ 



Ca I 

4227- 

9.3(0) Rayleigh 

l.l(-l)*cnf 3 

SW-ME 

“j" 

4044 

5.9(0) 

. 1. 7 (-1) 

SW-ME 

Ca II 

3934 

1*2(1) 

8. 3 (-2) 

SW-ME 

Fe I 

3441 

6*1(-1) 

1.6(0) 

SW-ME 

Ni I 

3370 

6.5(0) 

1.5(-1) 

SW-ME 

Ti I 

3342 

3.4(1) 

2. 9(-2) 

SW-ME 

***Na I 

3302 

2.6(0) 

3.9(-l) 

SW-ME 


3233 

3.8(0) 

2. 6 (-1) 

SW-ME 

kZ I 

3082 

8.7(0) 

l.l(-l) 

SW-ME 

Mg I 

2852 

1*2(1) 

8. 3(-2) 

SW-ME 

Mg II 

2796 

3.8(0) 

2* 6 (-1) 

SW-ME 

Fe II 

2599 

2.1(0) 

■4. 8(-l) 

SW-JIE 

Si I 

2514 

6.3(-l) 

1.6(0) 

SW-ME 



XX 500-2500$ 



C I 

1657 

1. 4(-2) 

7*1(1) 

SW-ME-VO-FC 

Xe I 

1470 

9.2(-5) 

1*1(4) 

SW-RA-PR 

• 0 I 

1304 

9. 0(-5) 

1.1(4) 

SW-ME-VO-FC 

’H I 

1216 

4. 9(-l) 

2.0(0) 

SW-VO-FC 

N I 

1200 

2. 6 (-4) 

3.8(3) 

SW-ME-FC 

Kr I 

1165 

4. 9(-6) 

2.0(5). 

SW-RA 

Ar I 

1048 

l.l(-5) 

9.1(4) 

. SW-RA-PR « 

Ne I 

736 

4.0(-6) 

2.5(5) 

sw 

He I 

584 

I. 8(-4) 

5.6(3) 

SW-RA 

^’number 

in parentheses denotes power of ten, i.e. 1.0(2) - 

* 1.0 x 10 2 



*’*SOURCE FUNCTIONS: 

SW, solar wind; ME, meteoric bombardment; VO, volcanic activity; RA, radioactive decay; 

PE; original atmospheric residua; FC, fuel contaminants 
the tabulated wavelengths pertain to second resonance lines of the species; the first resonance 
lines are K I (7665), Na I (5890). Li I (6708). 



dichotomy. The instrumentation is duscussed in greater detail in the 
following section where it is demonstrated that system detectabilities 
can be achieved which correspond to signal intensities between about 1 
to 10 Rayleighs. On this basis, the data in Table 14 indicate that this 
instrumentation capability will result in high probability of detecting 
and identifying (or alternatively, establishing meaningful upper limit 
surface density values) a large number and variety of possible lunar 
atmospheric species. 

The proposed experiment is predicated on the employment of a simple 
design spectrograph of reasonable resolution which is both compact and 
' light weight. An optimum design for the present experimental purpose is 
shewn in Figure 71 which illustrates the optical configuration of a simple 
normal incidence spectrograph. In Figure 71 the housing geometry is in- 
dicated in order to illustrate the compact nature of the resultant 
spectrograph. The entrance slit,’ grating and plate holder positions in- 
dicated in the figure lie on the Rowland circle so as to achieve suffic- 
iently good focus for an extended wavelength region. An existing proto- 
type for the present application has been designed, fabricated, and cal- 
ibrated by McPherson Instrument Co. (a division of GCA Corporation) and 
utilized in a variety of rocket experiments designed to obtain vacuum 
ultraviolet spectra in the earth upper atmosphere is shown in’ Figure 72. 

A technique of stacking the two module assemblies for compactness is also . 
illustrated in the figure. Further details of the basic spectrograph 
instrument are not included herein since these are avaife ble in the lit- 
erature (Ref. 132). The recommended optical configuration involves over- 
all unit assembly spectrograph dimensions of 4 inches depth by 8 inches 
width by 10 inches length (i.e., 320 in^) with a weight of less than 5 lbs./ 
instrument. Finally, it can be seen from the schematic of -Figure 72 that 
the option is available to incorporate either a single spectrograph or 
the dual instrument on any specific Apollo mission. 

The critical instrumental parameters are listed in Table 15 for both 
spectrograph designs wherein the indicated values have been assigned ten- 
tatively on the basis of a preliminary analysis of the present problem. 

It should be stressed at this point that all of the parametric require- 
ments stipulated in this tabulation can be accomodated by off-the-shelf 
and state-of-the-art components. Brief discussions are presented below * 
of the individual parameters included in Table 15 for both instrumental 
designs. 

The indicated gratings are available from Bausch and Lornb (Ref. 133) 
and other suppliers. The overall instrumental transmissivity estimate of 
Of 5 for XX 25 00-45 00A is considered conservative owing to the high reflec- 
tivity of materials and grating efficiencies of this spectral region. 
Alternatively, for the spectral region XX 'less than 1500&, it is well- 
established (Ref. 133) that reflectivities can be enhanced significantly 
by employing stock item platinum coated gratings which have been over- 
coated with nagnesium fluoride. Therefore, an overall transmissivity of 
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figure 71 An optical configuration of a normal incidence spectro- 
graph. The housing geometry is shown to illustrate the compact 
nature of this configuration. The dotted portions represent extension 
of Rowland circle. 






TABLE 15 

INSTRUMENTAL PARAMETERS FOR PORTABLE -SPECTROGRAPHS 


Spectrograph #1 
AX 2500-4500$ 


Spectrograph #2 
XX 500-250 q8 


Optical Specifications 


Grating Specifications 
Grating Type 

Focal Length, cm 
Grating blaze, $ 

Grooves/mm 
Ruled area, mm 
Reciprocal dispersion i?/mm 
Omega, ft, steradians (xl0~) 

Slit width, cm 
Slit length,, cm 
Overall Transmissivity 

Film Specifications 
Typ, e 

. Wavelength sensitivity .A 
Exposure Required for 
minimum density E, 
ergs/cm^ (x lCr) 2 

Mean E value, photons /cm 
. (x 10“ 7 ) 

Reciprocity Failure 
Availability 

Instrument Overall Specifications 

Size, inches 
Weight, lbs. 

Power requirements 
. Data measurements requirements 
Data support requirements 


Concave blazed 
reflecting 
20 

3200 (first order) 

1200 
53 x 56 
«-40 

7.4 (FOV ~15°) 

0. 10 
1.0 
0. 5 


Kodak 2485 

2200-6700 

5.0 


(X = 3200R) 
Negligible 
Roll Form 


4 x 8 x 10 

5 

None 

None 

None 


Concave blazed 
reflecting 
20 

1200 (first order) 
600 (second order \ 
1200 
53 x 56 
- 40 

7.4 (FOV ~15°) 

0 . 10 
1.0 
0.5 


Kodak j-ui-ujL 

50-4000 

4.2 


2. 5 

(x = 1200X) 

Negligible 
Roll Form 


4 -x 8 x 10 

5 

None 

None 

None 
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0.5 was selected throughout the entire spectral region for the purposes 
of the present discussion. In practice, it is planned to establish the 
transmissivities of the individual spectrographs on the basis of a lab- 
oratory measurement program. 

The rationale for the film selections of Table 15 can be discussed 
in terms of Figures 73, 74, and 75. Specifically for XX 2500-4500$ Kodak 
2485 film was selected on the basis of the spectral sensitivity curves 
of Figure 73 which was extracted directly from the appropriate Kodak data 
sheet. These published data were employed to obtain the corresponding 
E-values shown in Table 15. In the present experimental mode, it can be 
anticipated that relatively long exposure- times may be involved so that 
the role of reciprocity failure must be considered. In this regard, it 
should be stressed that Lewis and James (Ref. 134) have concluded re- 
cently that reciprocity failure was negligible when the film was employed 
under vacuum conditions. The resultant behavior is summarized in Figure 74 
which was extracted directly from the cited paper. In addition, the 
results of the figure indicate that an increase in sensitivity occurs 
under vacuum conditions although this conclusion should be considered 
tentative. In any case, for the present purposes reciprocity failure was 
assumed to be negligible so that the sensitivity curves of Figure 73 are 
applicable for exposures up to about 10^- seconds. 

For the spectral region XX 500-2500$, the Kodak 101-01 film sel- 
ection was predicated on the following rationale. An investigation was 
performed recently by W. K. Fowler, et al (Ref. 135) on the absolute cal- 
ibration of SC-5 film employed for far ultraviolet rocket spectroscopy. 

A series of exposures of SC-5 film at a wavelength of 1200$ is shown in 
Figure 75 which has been extracted directly from the cited paper (Ref. 135). 
It can be seen that, even for an exposure time of 1/20 of one- second and 
a beam intensity of 1.7 x 10“2 ergs/second, sufficient film exposure was 
obtained (for the indicated development process) as shown in the lower 
left slit image in Figure 75. On the basis of more recent experimental 
results, it is now recognized (Ref. 136) that Kodak 101-01 film is both 
about 50 percent faster than the SC-5 film and also results in a clearer 
.background -upon proper development. Additionally, it has the further 
advantage of being available in roll form. On the basis of this in- 
creased capability, as well as the degree of darkening observed in the 
lower left slit image of Figure 75 estimates were derived for the exposure 
. required for the minimum operating density of Kodak 101-01 film which .are 
presented in Table 15. With respect to reciprocity failure, Fowler, et al 
(Refs. 135,136) have concluded that for SC-5 and 101-01 film, it is com- 
pletely negligible. Therefore, it is appropriate to apply the tabulated 
E-value to the calculation of exposure times required for the purposes of 
the present experiment. 

It is suggested that the dual assembly spectrograph be operated within 
the protective shadow of the LEM. The optical axis of the instrument should 
be pointed about 10° above the lunar horizon with the slit oriented in 
the anti-solar position in order to minimize any degradation introduced 
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LOG SENSITIVITY— Reciprocal of Exposure (ergs/cm 2 ) 
Required to Produce the Densities Given 



FIGURE 73 

Spectral sensitivity of Kodak 
2485 high speed recording film. The 
dotted curve obtains for D =~ 0.10. 





0 0 1.0 2 0 

Log relative E 


fxcurs 74 characteristic curves 
for blue light exposure in room 
air and vacuum of pure bromide 
emulsion without deliberate 
chemical sensitization. Expos- 
ure times are noted on each curve 




Figure 75 - Series of exposures on SC- 
film at 1200&. Exposure tines 
were 1/20, 1/10. 1/5. 1/2- 1, 2, 
4. & 8 sec to a beam of 1.7 x 
10“ 2 ergs /sec total flux. After 
exposure the film was cut in two 
and the upper half developed in 
OK- 20 for 12 min while the lower 
half was developed in D - 1 9 b tor 
2 minutes- 
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by surface scattered radiation- This observational configuration can be 
instituted by the astronaut either by means of a built-in tripod arrange- 
ment or by simply suspending the dual spectrograph from a pre-selected 
location on the bottom section of the LEM vehicle * The astronaut wou Id 
then initiate the instrumental operation mode wherein the entire sequence 
of pre-programmed exposures would be performed automatically through 
equipment shutdown (which occurs after several hours of operation), The 
only other required astronaut activity involves removal of the portable 
roll rilm holder for return to earth for subsequent film processing and 
data reduction and analysis. Thus, it is evident that although minimum 
astronaut activity is involved in the conduct of the proposed experiment , 
his role is vital to its successful performance* 

As indicated previous ly, meaningful const! tuence data will be ob- 
tained for a system detectability of between 1 and 10 Rayleighs as shown 
in Table 14. On this basis then, a number of suggested exposure times * 
have been selected employing the parametric instrumentation values dis- 
cussed above and listed in Table 15* At this point, it is convenient to 
define the following useful parameters: 

^ - the brightness value for the background radiation con- 

tained within the spectral resolution of one slit 
element- For convenience, this is given in units, of 
Rayleighs (see Figure 76). 


X g = the brightness for the spectrally discrete signals falling 
on the entrance slit* For convenience this is given in 
units of Rayleighs (see Table 14} 

2 

A “ slit area, cm 


(1 pj viewing solid angle, steradiaus 

T - overall transmissivity, 0*50 

t - time of exposure, seconds 




E 


integrated exposure relevant to one slit area element image 
on photographic film, ergs /cm^- slit area and/or photons/ 
cm^-slit area 


't/A; integrated exposure relevant to 1 cm^ image of photo- 
graphic film (see Figures 75 and 74). 


It is appropriate now to evaluate the exposure times required to 
achieve a I to 10 Rayleigh detection capability. This can be estimated 
in the following manner : 
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(Xg+Xg) as a function of wavelength for 
mean galactic sky conditions and average 


QZS-ti’OtO 


[(X s + X B ) t) 


(53) 


$ = (An.) r 




so that 

ft T [X s + X B ) c ) 
E = 4n 


substitution from Table 15 yields: 


(54) 


E = 3.0 x 10" 3 [(X e + X ) c] (55) 

S D 

Equation (53) is plotted in Figure 77 where E-values of between 10^ and 
10^ photons/cm^ (see Table 15) are plotted as a function of exposure 
times (t) in the range between 10 and 300 minutes. For future reference 
a number of selected (Xs + X B ) -values are included which range between 
0.5 and 40- Rayleighs. On the basis of these data, the background data 
shown in Figure 76, and the instrumental characteristics listed in. Table 
15, it will be demonstrated that a 1-10 Rayleigh signal brightness de- 
tectability can be achieved by proper operation of the dual spectrograph 
on the lunar surface. 

The results of Figure 77 demonstrate the exposure time is selected 
on the basis of the anticipated (X g + X B ) -value. It is important to 
stress here that the X B -value is pertinent to the determination of the 
proper exposure time. Specifically, it is evident from the data in Fig- 
ure 76 that the Xjj-value exhibits significant variation over the spectral 
regions suggested' for the two spectrograph instruments. Specifically, 
for XX 2500-4500$, X B -values o„f between 4 and 40 Rayleigh (spectral 
resolution ~ 40$) obtain. Alternatively for the spectral region XX 500- 
2500$, the background contribution can be ignored for the most part. 

For this reason, exposure estimates are obtained individually for the 
two spectrograph systems. 

With respect to XX 2500-4500$, a practicing rule of thumb of photo- 
graphy involves successful image detection against a strong background 
when the brightness of the image amounts to at least 10 percent of the 
background, i.e. (X s + X B )/X B > 1.1. However, this rule prevails for 
photographic systems operating on the straight portion of the film gamma 
curve. This is illustrated in Figure 78 where X s = X B for which the 
resultant discrimination of signal against background is readily apparent. 
Owing to the lack of relevant information on the gamma -curves of the pre- 
sently considered film systems, a conservative approach is adopted wherein 
it is assumed that spectral discrimination can be achieved when (X s + X B )/ 

X B > 1.25. On this basis, it is evident that proper exposure time selection 
must be predicated on the basis of the magnitude of the applicable back- 
ground brightness rather than the signal brightness for XX 2500-4500$. 

Thus, according to Figure 77, for an E-value of 8.0.x 10' photon/cm 
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-VALUES (XIO"') photons /cm 


0C239-I0I 



t -values (min) 


2 

figure 77 Plot of E-values in photons/cm vs_. t-values in minutes 
for a number of selected (X + X ) - values. These data can be 
employed to relate r<?quired S exposure times to desired detectivity. 



figure 78 Schematic 
background film 


of typical photographic 
densities are equal. 


spectra wherein the signal 


and 



presented in Table 15, exposure times of between 10 and 100 minutes are 
required for Xg-values of 40-4 Rayleighs, respectively. According to 
the detection criterion discussed above, a corresponding signal detection 
capability would be realized of between 10 and 1 Rayleigh, respectively. 

As a consequence, a recommended program of exposures would be employed 
in practice which includes exposure times of 1-300 minutes. 

In contrast to the above situation, the spectral region XX 500-2500$ 
involve no .significant background contributions according to Figure 76. 

For this case, exposure time selection is predicated only on the basis of 
the magnitudes of the' expected X s -values. In this regard, the data of 
Figure 77 suggest that for the derived E-value of 2.5 x 10 7 photons/cm 2 
similar exposure times as those derived above for the XX 2500-4500$ spec- 
tral region obtain for a 1-10 Rayleigh signal detection capability for 
the XX 500-2500A spectral region. Accordingly, it is recommended that an 
exposure time sequence identical to that obtained previously for XX 2500- 
- 4500A of between 10 and 300 minutes be selected for this lower spectral 
region. 

Relevant to the above discussion, it. may be noted that a well-estab- 
lished practice involves the controlled pre-exposure (to threshold) of - 
fast film in order to enhance its speed characteristics. 

The above suggested exposure time sequences should be regarded as 
tentative estimates since the final selection will be determined both on 
the basis of the specific prevalent background conditions as well as 
laboratory calibration of the optical components and films in the indiv- 
idual spectrograph systems. In this regard, an alternative instrumental 
configuration which would involve the employment of image intensification 
has not been considered herein owing to the additional instrumental com- 
plexities which would be introduced. However, in the event that such 
lower exposure times are necessary, this alternative could be applied to 
this experiment. 

The photographic spectra obta inecl in the conduct of the proposed - 
experiment require minimum analysis, since the film format intrinsically- 
yields the unambiguous spectral locations of resonance scattering of the 
individual species. Alternatively, detailed interpretation of these 
results is expected to involve a considerable difort. For example, it is 
anticipated that the explanation of positive results will involve a de- 
tailed investigation to identify the major processes contributing to the 
existence of specific species in the lunar atmosphere. In addition, any 
measured time variation of the content of a given species will involve 
the identification of specific phenomenological sink and source functions 
such as solar wind activity, meteoric impact, LEM fuel contaminants, etc. 

It should be emphasized, however, that even negative results can be employed 
to establish meaningful upper limit values for a large number and variety 
of possible lunar atmospheric constituents. These data can then be utilized 


•160 



to preclude or evaluate the roles of several suggested mechanisms con- 
tributing to the lunar atmosphere. Finally, the overall experimental 
results can be employed to evaluate the validity and feasibility of estab- 
lishing a permanent astronomical -air glow lunar surface station to monitor 
selected lunar atmospheric species and communicate these data routinely 
to earth. 

Successful performance of the proposed experiment would be assisted 
greatly by the performance of parallel, theoretical and laboratory inves- 
tigations. The theoretical effort should involve an investigation of the 
relationship between specific lunar atmospheric species and appropriate 
source and sink functions as well as the mechanisms involved in generating 
steady-state concentrations. The laboratory investigation should involve 
the measurement of the optical characteristics of the two suggested spec- 
trograph systems as well as absolute calibration of the films for employ- 
ment herein. These supporting data would then be employed in the cogent 
interpretation of the derived experimental results. 
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. III. MISCELLANEOUS TOPICS 


This section contains brief summaries of the scientific material 
generated under the present program and presented at scientific and/or 
professional meetings. In addition, there is included other miscellaneous 
topics of interest in the performance of the current contract commitments. ' 

- During the current reporting period, Dr. James A. R. Samson attended 
the IAGA symposium on "Laboratory Measurements of Aeronomic Interest" 
held at the University of York, Toronto during 3-4 September 1968 and also 
to visit the laboratories of Dr. Marmet at the University of Laval, Quebec 
■ City, Canada. 

Of particular interest to the present contract was the paper presented 
by Dr. R. Schoen on Photoelectron Spectroscopy, Photoionization and Photo- 
absorption Measurements. He pointed out that the absorption cross sec- 
tions of CO 2 vary widely from observer to observer and although he believed 
the values obtained at GCA were the best available he felt that they should 
be remeasured. He also brought up the importance of the angular distribu- 
tion of the photoelectrons in photoelectron spectroscopy experiments. 
Independently of the work produced at GCA, he also arrived at the look 
angle of 54 44* for observations in order to eliminate any confusion of the 
results due to varying angular distributions. 

McElroy brought out the fact that the abundance of atomic nitrogen 
produced by dissociative ionization appears to be too small to account 
for what he expects to exist in the upper atmosphere. The discussions on 
this point suggested that electron impact collisions with N 2 may be the 
major instrument for producing dissociative ionization. 

Dr. Hoffman of the AFCRL, gave an excellent review of the present 
state of - absorption cross sections of ‘atmospheric gases and pointed out 
the need for continued measurements at higher wavelength resolution. 

At Laval University, Dr. Marmet described his new ideas for measuring 
the o kineti c energies of photoelectrons. His new analyzer is basically a , 
127 cylindrical deflection analyzer with a modification to provide a 
higher transmission of electrons. The device looks extremely promising 
for photoelectron work. However, it has not been tested as yet and we 
propose to hold off making this type until such time as Dr. Marmet has 
more information on its energy resolution and other characteristics. Dr. 
Marmet has had many years of experience dealing with electron and ion op- 
tics and he was able to provide invaluable suggestions about the techniques 
in working with electrons, especially in methods to remove local magnetic 
fields. 

Dr. James A. R. Samson attended the 21st Annual Gaseous Electronic 
Conference held at the University of Boulder, Boulder, Colorado, 16-18 
October 1968. Dr. Samson presented a paper entitled "Higher Ionization 
Potentials of Molecules Determined by Photoelectron Spectroscopy. " This 
paper was received with considerable interest since there is an accelerating 
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activity in the new field of electron spectroscopy. Historically, this 
conference deals with photon -atom in terse tion, a tomcat cm inters c cions, ion- 
molecule reactions and electron-molecule interactions. Many of these 
areas are of direct interest to the performance of the current contract 
commitments. Of particular interest to some current problems in photo- 
electron spectroscopy was a paper entitled "Configuration Mixing of 
Continuum States* 1 presented by A. J. Mendez. This theoretical paper has 
been of significant importance in defining future pertinent experimental 
measurements. For example, it x^as" pointed out that it would be of sig- 
nificant importance to measure the coefficients A and B for the angular 
distribution of photoelectrons in wavelength regions wherein auto- ioniza- 
tion, and overlapping continuum states' occur- With the experimental values 
availalbe, the theoretician could then determine the proper wave functions 
of the states ■ involved. 

During the conference local visits were made to NCAR and E$SA* E. 
Ferguson (E5SA) deomonstrated laboratory evidence for the reaction He + + 

N 2 ** .N2"^ ■h.hv (Rf.2000^). This was of special interest to our laboratory 
effort, since GCA has previously observed "fluorescence" in He at wave- 
lengths shorter than those corresponding to the ionization potential- On 
the basis of Ferguson *s data it would appear then, that these emissions 
could have been due to N2 impurities in the gas sample. 

Dr- Frederick F. Memo attended the GOSPAR meeting at Prague, 
Czechoslovakia, for the purposes of presenting a paper entitled, "On the 
Presence of Carbon Atoms in the tipper Atmosphere of Venus." A rather large 
number of other papers presented at the meeting were relevant to the task 
requirements included in the present contract. As such a significant 
amount of valuable input data were obtained and will be' incorporated in 
forthcoming reports as required. 


The the First Annual Meeting of the Division of Electron and Atomic 
Physics of the American Physical Society on 17-19 November 1969 a paper 
entitled "Angular Distributions of Photoelectrons and Partial Photo ionisa- 
tion Cross Sections" was given by Dr. James A. R. Samson of GCA Technology 
Division, Bedford, Massachusetts. An experimental arrangement was given 
for measuring the angular distribution of photoelectrons. The asymmetry 


coefficient, p, was found to be 0..3 for argon, and 0.3 and 0.12 for molec- 
ular nitrogen at 584^ when the residual ion core is left in its X £ and 
A 2 H states, respectively. The partial photoionisation cross sections were 
given for t>£ at 186, 209, 247, 330 and 46G&, New ionisation potentials 


for O 2 were found at 23.5, 24. 6‘ and 27.3 eV j- 0.3 eV. 
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